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M335 5.0/5.1 
Introduction 


Any function whose domain is the set of positive integers is called a number- 
theoretic or arithmetic function. Most of the number-theoretic functions which 
we shall meet also have the set of positive integers as codomain. 


А number-theoretic function f which has the property that f (mn) = f (m) f (n), 
whenever m and n are relatively prime integers, is said to be a multiplicative 
function. It is a simple inductive consequence of this definition that if f is 
multiplicative and n has prime factorization n = p*'p5? ... p* then 

f(n) = (pi) f (p?) ... f(p). Thus, to determine a multiplicative function 
completely, we require to know its value at each prime power. 


There are numerous interesting multiplicative functions which arise from 
looking at the divisors of an integer n. For instance, for any non-negative 
integer i, the function | 

fin) = ў 4, 


where the summation is taken over all the positive divisors d of n, сап be shown 
to be multiplicative. The first two of these, fo and fı, are very famous functions. 
fo, Which we rename т, is the function which counts the number of positive 
divisors of each positive integer n > 1. fı, which we rename o, is the function 
which adds the positive divisors of each integer n = 1. These two functions t 
and c and their properties are discussed in Section 5.1. 


Section 5.2 investigates properties of the greatest integer function. Although this 
function is not number-theoretic, because of its importance in number theory it 
holds a natural place in this unit. 


Probably the best known of all multiplicative functions is Euler's $-function, 
which assigns to each positive integer n the number of integers in a complete set 
of residues which are relatively prime to n. Section 5.3 introduces the $-function 
and proves that it is multiplicative. Section 5.4 looks at properties of the ф- 
function and presents Euler's generalization of Fermat's Little Theorem, a result 
which is arguably one of the most useful in the whole of number theory. 


The student with some background in algebra may well recognize that much of 
the material in Sections 5.3 and 5.4 could be set in the wider context of group 
theory. We have made several remarks to this effect in the unit. However, the 
emphasis of the course remains with ‘numbers’ and these algebraic asides can be 
ignored safely. 


The unit comprises Sections 6.1, 6.3, 7.2 and 7.3 of B. 


The Functions t and o 


In this section we introduce the idea of multiplicative functions and meet our 
first two examples: the functions т and о. 


Theorem 6-1 of B, in the first reading passage, ought to be familiar, as this 
result occurred in Unit 2 when we were looking at the canonical prime 
factorization of positive integers. 


READ B Section 6.1, page 110 as far as page 114, line — 7. 


(i) Page 110, line —5 
By definition, an integer n > 1 is prime if and only if it has exactly two 
positive divisors, namely itself and 1. That is, n is prime if and only if 
t(n) = 2, or alternatively о(п) = n + 1. 
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(ii) Page 112, line —4 


When the product is expanded, the result is a sum of terms each of which 
is a product of r factors, one from each of the brackets. (In fact, there will 
be (К, + 1) (к +1)...(k, + 1) = t(n) terms.) In other words, the product is 
the sum of the t(n) distinct terms of the form рр: - + p¢, where 

0 < a; € k; and by Theorem 6-1, this sum is precisely o (n). 


(iii) Page 114, line 14 


For an illustration, take n = 18. As d runs through the т(18) divisors 1, 2, 
18 

3, 6, 9 and 18, so d' — p takes the values 18, 9, 6, 3, 2 and 1 respectively. 
а 

So i 
Па= Па =1-2-3-6-9-18, 


d|18 d'|i8 


Multiplying the т(18) = 6 equations of the form 18 = dd’ together gives 


18° = I Ш) = (1:2:3:6:9:18)2 


alts] (ав 


showing that the product of the positive divisors of 18 is equal to 18°. А 


As an example, let us look at the Problem quoted. 


EXAMPLE 1 B page 118, Problems 6.1, number 7. 


(a) We use Theorem 6-2 and SAQ 4 of Unit 2. If n is a perfect square, then the 


(b 


— 


exponents in the prime factorization of n are even, i.e. n = p?" p^ ... p2, 
Therefore 
a(n) = (2l, + 1) Ql; + 1): I, +1) 


is a product of odd integers, which is odd. 


On the other hand, if n is not a perfect square, then at least one exponent in 
the prime factorization n = pp? +-+ p is odd, so at least one factor in the 
product 

(ky + 1) (ko + 1): (k, + 1) = (n) 
is even, and hence т(п) is even. 
(Alternatively, without appealing to Theorem 6-2, we could observe that if d 
А a ae n Р - 
is a divisor of n then so is 7 since n = d: i} This means that the divisors 

а 

of n fall naturally into pairs, and hence т(п) is even, unless there is a divisor d 


T pge i ; ? 
of п such that d = Т This in turn means п = d", i.e. n is a perfect square.) 


Let n have prime factorization n = pf'p? ---p*. The formula at the foot of 
page 112 of B gives 
(п) = [| (1 +p, +p? ++ рр). 


1<i<r 


For a(n) to be odd, each factor in this product must be odd. Following the 
given hint, if p, is odd then 1 + p; + p? +--+ + pf is odd precisely when К, is 
even, whilst if p; = 2 then this term is odd for any value of Ку. Therefore c (n) 
is odd if and only if the exponent of every odd prime in the factorization of 
n is even, that is, n = 2*p2! p25 - -- p?' But then, if k is even, 


2۴2 2 


n = QV? pps pir), a perfect square; 


whilst if k is odd, 
kmt 
n=2(2 2 pep --- pr), twice a perfect square. 


Hence oc (n) is odd precisely when n is a perfect square or twice a perfect 
square. 


SAQ 1 
SAQ 2 


SAQ 3 
SAQ 4 


SOLUTION 
TO SAQ 1 


SOLUTION 
TO SAQ 2 


SOLUTION 
TO SAQ 3 


SOLUTION 
TO SAQ 4 


NOTES 
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B page 118, Problems 6.1, number 8. 


B page 119, Problems 6.1, number 9. 
(n is ‘square-free’ if it is not divisible by the square of any integer greater than 


1) 
В page 119, Problems 6.1, number 12. 


If t(n) = p?, where p is prime, describe the possible prime factorizations of n. 
Hence find the smallest positive integer n for which (i) t(n) = 8 and 
(ii) t(n) = 27. 


Let the positive divisors of n be dı, d>,...,d,. Then 


йй, 1 
==> tate te 
2 А d d, 
JE gt Seng 

_& dy d, 


n 


a(n) 


n 


(As d runs through the divisors of n, so does > and so the numerator is 
o(n).) E 


If n is square-free it is not divisible by the square of any prime, so the exponents 
in the prime factorization of n are all 1, ie. n = pp; ^: p, for distinct primes 
Di P2- --, D. Theorem 6-2 (a) now gives t(n) = 2" as required. В 


(a) If n has prime factorization n = pp? --- pk, then Theorem 6-2 (a) shows 
р 


that t(n) is the product of r factors of the type (k; + 1), each of which must 
be at least 2. As 10 cannot be expressed as the product of more than two 
such factors, we know that t(n) = 10 implies that r = 1 or r = 2. 

When r = 1, t(n) = Кү + 1 = 10 gives К, = 9. So any integer п of the form 

п = p? has т(п) = 10. 

When r = 2, t(n) = (Кү + 1) (k +1) = 10 = 2-5 gives Кү = 1 and К, = 4 (or 
vice-versa). So any integer n of the form n = p,p$ has t(n) = 10. 

The smallest value occurs when n = 3-2* = 48. 

(b) Using the hint, we know that if c (n) = 10 then 1 < n < 9, so it remains to 
check only that c (n) ¥ 10 for these nine values of n. The successive values of 
o(n) are 1, 3, 4, 7, 6, 12, 8, 15 and 13. BM 

There are three ways in which p? can be expressed as a product of factors each 


greater than 1, namely p: p: p, p: p? and p°. Accordingly, t(n) = p? arises when n 

is one of the forms pf "pp !, pf ! pZ^^! and p?! (this by Theorem 6-2 

(a)), where pi, P2, рз are distinct primes. 

G) t(n) = 8 = 2? arises when n = pipaps, or n = p,p} or n = pj. The least 
integer n of one of these forms is 3:2? = 24. 

(й) c(n) = 27 = 3? arises when n = pipipi, or n = pip$ or n = p25. The least 
integer п of one of these forms is 22-32-52 = 900. Ш 


READ B page 114, line —6 to the end of Section 6.1. 


(i) Page 115, line 12 
The use of f in B might be slightly confusing here. On the previous line, he 
refers to the specific function given by f(n) = 1 for all n > 1, observing it to 
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be a simple example of a multiplicative function. Now we revert to an 
arbitrary multiplicative function f, showing how Definition 6-2 can be 
extended to the product of any number of positive integers which are 
pairwise relatively prime (that is, for which any two are relatively prime). 
(ii) Page 116, line 15 

In the statement of the Lemma, B includes the hypothesis ‘and 

ged (d,,d,) = Г. Whilst it is important to appreciate that the divisors d, 
and d, are relatively prime, the inclusion of this statement is unnecessary 


since it follows automatically from the hypotheses ged (m,n) = 1, d,|n and 
dm А 


EXAMPLE 2 B page 120, Problems 6.1, number 22. 


(a) Putting s = 0, and 1, respectively, we have 


со(п) = Yd? = Y 1 = т(п), 
an 


d|n 
and 
e,(n) = Yd! = Уа = о(п). 
d|n d|n 


(b) The function f defined by f (n) = п? is multiplicative since 
J (mn) = (тп) = mi? = f (m) f (n) for all integers m and n, which certainly 
includes the case of m and n being relatively prime. (Functions like this 
which possess the multiplicative property for all pairs of positive integers are 
called completely multiplicative functions.) 


We can now apply Theorem 6-4 to see that 


F(n-»f(d)- Dd = о,(п) 
d\n d\n 
is multiplicative. 


(c) Let p* be any prime power. The positive divisors of p* are 1, p, p?,..., p* and 
so 


(РУ) = 18 + р? + p+ +++ +p 


E (p+: =Ñ 
=i 


Using the multiplicative property, 


(by SAQ 2 of Unit 1, with a = 1, r = p’). 


a(n) = as(pi') o(p?) o,(pF) 


a = ') = MES ') A (= m ) 

pi-t 5 —1 peat 
As a numerical illustration of this formula, let us determine the sum of the 
cubes of the divisors of 40. 


°з@)- 902-8) = (= S 
= (2° + 1)-(23 + 1): (53 +1) 
(using a?" — 1 = (a" — 1) (a” + 1) three times) 
= 73710. 


SAQ 5 B page 119, Problems 6.1, number 14. 
SAQ 6 B page 120, Problems 6.1, number 20. 


SOLUTION 
TO SAQ 5 


SOLUTION 
TO SAQ 6 
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(a) From Theorem 6-2 (b), if n = 2*^! 
-D+ 1 


a(n) z =7—1=2л-1. 


(b) If 2* — 1 is prime then o(2* — 1) = (2* — 1) + 1 = 2*, and from part (a), 
c (2^!) = 2* — 1. Now since рса (2^ !, 2* — 1) = 1, 


a(n) = 6)2۴ !)eQ* — 1) = Q* — 1)2* = 2n. 
(An alternative method of proof would be to use Theorem 6-2 (b).) 
(c) o(n) = oQ* ^ 1)0 Q* — 3) (since ged (2*71, 2* — 3) = 1) 
= (2* — 1)(2* - 3 + 1) (since 2* — 3 is prime) 
= 2k(2* —3)-222n42. W 


(a) First note that the definition of р may be rephrased as p(n) = 2” where r is 
the number of distinct primes that divide n. Now if gcd (m,n) = 1, then there 
is no prime which divides both т and п, so if r distinct primes divide п, and 
s distinct primes divide m, there are r + s distinct primes dividing mn. 
Therefore 

p(mn) = 27** = 272" = p(m) p(n), 


showing p to be multiplicative. 


(b) F(n EU p (d) is multiplicative, by part (a) and Theorem 6-4. Now let n 


have sind factorization n = pf'p? +++ p*, so that 


F(n) = F(pt) F (p2) F(pr^) 


x | z | у, p) - it > ғи). 


dlp dip 


Our remaining task is to calculate 2;^ p(d) for a general prime power р“. 
d|p* 
Now the only divisors of р“ are 1, p, p?,..., p* so 


Уу p(d) = p(1) + p(p) + p(p?) +++ + p(p") 


d|pk 


-1-2-42-4:-5-2 
(since p is the only prime dividing p', r > 1) 


= 2k +1. 
Therefore 


F(n) = (2k; -1)Qk,—1)-Qk + 1). 8 


In SAQ 5 (b) we found that all integers n of the form n = 2*~1(2* — 1), where 
2* — 1 is prime, have the property that o(n) = 2n. That is п(=с(п) — n) is equal 
to the sum of all its positive divisors except itself; for example с(6) —12 or, 
alternatively, 6 = 1 + 2 + 3. Integers n with this property, с(п) = 2n, are called 
perfect numbers. The four smallest perfect numbers are 6, 28, 496 and 8128, 
which are all integers of the above form, 2* ^! (2* — 1), where 2* — 1 is prime, 
corresponding respectively to k — 2, 3, 5 and 7. 


Perfect numbers first aroused the interest of the Pythagoreans, and the result of 
SAQ 5(b) appears as Proposition 36 of Book IX of Euclid's Elements. The 
history of the search for perfect numbers, and the many conjectures about them, 
makes fascinating reading. There are no known odd perfect numbers and though 
considerable inroads have now been made into proving that there are none, a 
complete proof is still awaited. On the other hand, it was shown by Euler that 
the perfect numbers of the form 2*^ ! (2* — 1), where 2* — 1 is prime, are the only 
even perfect numbers (see SAQ 7). 


SAQ 7 


SOLUTION 
TO SAQ 7 
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This result is not the end of the matter, for we are left with the problem of 
deciding which values of k make 2* — 1 prime. The numbers M 00 = I are 
called Mersenne numbers and if M, is prime it is called a Mersenne prime. If k is 
composite, the Mersenne number M, is easily shown to be composite, using the 
fact that 2" — 1 divides 2" — 1, so the only candidates for Mersenne primes are 
the numbers M,, where p is prime. As commented above, for the first four 
primes, M, turns out to be prime also. However, the hypothesis that M p might 
always be prime is refuted by the fact that M,, is composite. Because of 
continued interest, over the years, in the problem of finding Mersenne primes, 
the largest known prime has always been a Mersenne prime. As of mid-1979, 
the largest known prime is M44497, which has 13395 digits. (Although it is 
known that there are infinitely many primes it is often difficult to decide 
whether a given large integer is prime or not!) 


Chapter 10 of B is concerned with the search for perfect numbers and Mersenne 
primes and, although we have omitted it from this course, if you have time we 
would recommend you to read it. The chapter is much concerned with the 
historical development and, omitting some of the technical details of proofs, it is 
very readable. In any case, to appreciate the difficulty of finding Mersenne 
primes, you might find the anecdote in the second paragraph on page 226 
stimulating. (Incidentally, Chapter 10 rounds off by mentioning another 
historically interesting family of integers: the Fermat numbers F „ = 27" +1. As 
Fı, F5, F5 and F, are all primes (and Е; is very large), the conjecture was 
voiced for some time that all Fermat numbers are prime. But to this day no 
other Fermat primes have been found and there is a suspicion that all 
subsequent Fermat numbers are composite—a glance at page 236, line —9, will 
illustrate the difficulty of finding factorizations of ‘small’ Fermat numbers.) 


(In this exercise, we shall prove that every even perfect number n is of the form 
n = 2*71 Q* — 1), where 2* — 1 is prime.) 


Let n be an even perfect number, and write n — 26-1, where k > 2 and m is 
odd. Show that the requirement that п is perfect (a(n) = 2n) gives 


2*m = (2* — 1)e (m). 


Deduce from this that (2* — 1)|т, and, putting M = , establish that M — 1 


MCN 
2-1 


and 2* — 1 = m is prime. 


See Theorem 10-1 on page 220 of B. The proof begins with the converse result, 
which is our SAQ 5(b). B 


The Greatest Integer Function 


Though not a number-theoretic function, the greatest integer function will prove 
extremely useful to us in Unit 6. For this reason, we follow B in studying this 
function next. We shall, however, omit those results for which we do not have 
immediate use. 


There is a slight conflict in that the notation used by B for the greatest integer 
function is the same as our earlier usage for congruence classes. The square 
bracket notation is quite standard for both these contexts. No confusion should 
arise, since congruence classes appear only when we are dealing exclusively with 
integers, whilst the domain of the greatest integer function is the set of all real 
numbers. 


[READ B Section 6.3, page 126 as far as the foot of page 128. 


NOTES 


SAQ 8 
SAQ 9 


SOLUTION 
TO SAQ 8 
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(i) 


(ii) 


(iii) 


(iv) 


(a) 


(b) 


Page 126, line 13 

In fact, [x] is the unique integer n which satisfies x = n + 0 with 

0 < 0 < 1, since this implies that x — 1 < n < x. 

Page 126, line — 4 

To illustrate the proof of Theorem 6-9, look at the particular case when 

n = 37 and p = 3. Our objective is to find the exponent of the highest 
power of 3 which divides 37! = 1۰2۰3۰ ‘°° -36-37. Since 3 is prime, we 
can obtain this power by adding together the powers of 3 which divide the 
factors 1, 2, 3,...,37. Of these factors, twelve are divisible by 3: 


3(=1:3), 6(=2-3), .., 33(=11:3), 36(=12:3). 


That there are twelve factors divisible by 3 comes from the fact that 
36 = 12-3 is the largest multiple of 3 which is less than or equal to 37; in 
other words, because 12 = [37 ]. 


So we know 3!2|37! However, four of the factors of 37! are divisible by 32: 
9(21:32, 18(=2:32), 27(=3:32), 36(=4:32). 
That there аге four such factors comes from the fact that 36 = 4-3? is the 


largest multiple of 3 which is less than or equal to 37; that is, 4 = [37]. 


Each of the above four factors contributes one more to the exponent of 
3 in 37! so that we now know 3!?*^/37! Finally, 27 = 1-3? is actually 
divisible by 3° and so contributes a further 1 = [37] to the exponent of 3 
in 37! So the total exponent is 12 + 4 + 1 = [2] + [32] + [53]. 


Page 128, line 7 
That [a + b] > [a] + [b] is one of the results we ask you to prove in SAQ 
8, which is intended to familiarize you with the greatest integer function. 


Page 128, line 11 
To illustrate formula (1), let us look a1 the example п = 37, р = 3, г = 11. 
The inequalities for the successive values of k are 


к= 1 RS [0] + [26]; 122348 
k=2 Bz [5] + [20]; 42142 
k=3 > 2474+ [28]; 1>0+0 


and for subsequent values of k every term is 0. Adding: 


п> у Bel + + 3 By 1724+10. 
k= k=1 


= 
Ш 
- 


Now, by Theorem 6-9, y [34] is the exponent of the highest power of 3 
к=1 


which divides 37! Similarly, the terms on the right-hand side of the 
inequality are the exponents of the highest power of 3 dividing 11! and 26! 
respectively. The right-hand side therefore gives the exponent of the highest 
power of 3 dividing the product (11!) (26!), namely 3**1°, and the 
inequality tells us that this is not greater than the highest power of 3 
dividing 37!, namely 3". А 


B pages 130 and 131, Problems 6.3, number 2 parts (a), (b), (c), (d) and (e). 
B page 131, Problems 6.3, number 5. 


Let x = [x] + 0, where 0 x 0 < 1: It follows that x + n = ([x] + n) + 0, 
where [x] + n is an integer and 0 < 0 < 1, and so [x + n] = [x] + n. 


If x is an integer, then so is —x, which means that [x] = x and 
[—х] = —х. Therefore [x] + [-x] 2 x — x = 0. If x is not an integer 


SOLUTION 
TO SAQ 9 


5.3 


M335 5.2/5.3 


then x = [x] + 0, where 0 < 0 < 1, and so, using the hint, 


x [x] —1+ (1 0). Now — [x] — is an integer and 
0«1—0 «1, therefore [Cx] = — [x] — 1, giving [x] + [-x] 2 — f. 
(c) [x] € x and [y] € y imply that [x] + [y] <x + y. Now [x + y] is the 
greatest integer less than or equal to x 4- y, and, as [x] + [v] is an integer 
and [x] + [y] € x + у, we have [x] + [y] € [x * y]. 
Writing x = [x] + 0 and y = [y] + ф, where 0 < 0, <1, 
xy = [x] [v]  0[y] + ¢ [x] + 09. Since x and y are positive, [x] > 0 and 
[y] = 0, and consequently 0[y] + $ [x] + 0ф > 0. Therefore xy > [x] [y], 
but since [xy] is the greatest integer which is less than or equal to xy, 
[xy] = [x] [y]. 
(d) Using the hint, we write = = H + 0 where 0 < 0 < 1, so that 


x=n Н + nO and hence by part (a) (since || is an integer), 


0 
[x] = “|= | + [n0]. Now 0 € n0 < п, so 0 < [n0] < n and 0 < BA PER 


х 5 ) i 0 
Hence Ta = Н + el: where Н is an integer and 0 < De «1; 
n n n n n 


therefore = = Н 
п п 


т пт 


Q 
= 
БЫ 3 
لا‎ 
lA 


an and hence n & = oo . E 
k k k 


Following Example 6-2, we could use Theorem 6-9 to find the highest power of 


each of the primes 2 and 5 dividing 1000! and the required answer will be the 
smaller of these two exponents. However, it is abundantly clear from the 
formula that 5 has the smaller exponent and so the number of zeroes in which 
1000! terminates is 

[22°] + [1929] + [02] + [4999] = 200 + 40 +8 +1 = 249. 
(If you care to check, you will find that the exponent of 2 is 994) Bi 


Euler’s $-function 


The Euler $-function is probably the most important of all multiplicative 


m m |. ; 
p therefore Н <——,80 2 is an integer less than or equal to 
К К К К 


functions. For each positive integer п, it counts the number of integers from 1 to 
n inclusive which are relatively prime to n. So, for example, ф(4) = 2 since 1 and 


3 are relatively prime to 4, and 4(5) = 4 since 1, 2, 3 and 4 are relatively prime 
to 5. 


COMMENT From the point of view of group theory, the $-function arises as 
follows. We have already seen that Z, is a group with binary operation +, so it 
seems natural to investigate Z, with binary operation x. We quickly see that a 
group does not result because, although x is closed and associative on Z, and 
[1] [x] = [x] shows that [1] is a multiplicative identity, not all elements have 
multiplicative inverses. For instance, since [0] [x] = [0] for all x, [0] does not 
have an inverse (except in Z, where [0] = [1]). In fact, the situation is worse. 
For example, іп Z4, [2] does not have an inverse since [2] [0] = [0], 

[2] 1] = [2], [2] [2] = [0] and [2] [3] = [2]. 


NOTES 
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But what about the subset U, of Z, consisting of those elements which do have 
multiplicative inverses, the so-called units of Z,? In this case, it is not difficult to 
show that the product of two units is a unit (closure), and then the associativity, 
identity and inverse axioms follow immediately. So U, is a group under binary 

operation x. 


Which elements are the units of Z,? [x] is a unit of Z, if and only if there is an 
integer a such that [x] [a] = [1], or what amounts to the same thing, 

xa = 1 (mod n). This, in turn, can be re-expressed as xa + nb = 1, for some 
integer b. Now Theorem 2-4 on page 27 of B shows that this is the case if and 
only if x is relatively prime to п. So there are ф(п) units in Z,, corresponding to 
the least positive residues of n which are relatively prime to n, and ф(п) is the 
order of the multiplicative group of units, U,,. 


READ В Sections 7.1 and 7.2, pages 134—140. 


(i) Page 138, line —6 
The proof of Theorem 7-2 may be clearer if we look at a particular case. 
If we take m = 9 and n = 8 the array is: 


(D 2 


e 
"iw 66 SSN 
£92605Q!zQo 
BAX RR RR a 
© 
OST OE OF © 
RALAR NAR 


64 65) 


In each column the entries аге of the form r, 9 + r, 2:9 n... 7:9 + n 
and are therefore congruent to each other modulo m — 9. So as far as 
being relatively prime to 9 is concerned, all the elements in a column 
behave the same way. We can therefore discard all the elements in the 
columns headed 3, 6, and 9, since those are not relatively prime to 9, and 
concentrate оп the ф(т) = $(9) = 6 columns headed by the integers 
relatively prime to 9, i.e. 1, 2, 4, 5, 7 and 8. 

From these remaining ф(т) columns, we now want to find which elements 
are relatively prime to n — 8. The key is in observing that, in each column, 
the n elements are congruent modulo п, to, in some order, 0, 1, 

2,...,7 — n — 1, since no two elements in a column can be congruent to 
one another modulo n. For example, in the column headed 4 we have 

4 = 4, 13 = 5, 22 = 6, 31 = 7, 40 = 0, 49 = 1, 58 = 2 and 67 = 3 

(modulo 8). It follows that in each of these Ф(т) columns, there must be 
(n) = ф(8) = 4 elements relatively prime to п, and these are the circled 
elements in the array. 


It follows that there are ¢ (m) b(n) = 6 x 4 = 24 elements relatively prime 
to mn — 72. 


(ü) Page 139, line 16 
Although B includes an inductive proof of Theorem 7-3, the result is 
immediate. From Theorem 7-2 we have 


$(n) = é(ptp? pF) = (pt) op): (pF) 


and Theorem 7-1 now completes the argument. A 


70 


EXAMPLE 3 


EXAMPLE 4 
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A word of warning. It is all too easy to remember 
(i) ф(р)=р—1, p prime 
and (ii) ф is multiplicative, 
and yet make erroneous deductions such as 
ф(р?) = ó(p): ó(p) = (p — 1). 


Make sure you remember the ‘pairwise relatively prime’ condition in the 
definition of multiplicative functions. 


To illustrate Euler’s ¢-function, we look at three examples. 


B page 140, Problems 7.2, number 4. 
(a) As n is odd, n and 2 are relatively prime. So 
$Qn) = ф(2)ф(п) = p(n), 

since ф(2) = 1. 

(b) Let n = 2'm where r > 1 and m is odd. 
$n) = ф(2'*1т) = $(27*') фт) = 2p (m), 
and 
26(n) = 26(2'm) = 262") $(m) = 2۰27 1p (m) = 2p (m). 

Therefore $ (2n) = 26 (n). 

(c) If 3|n, put n = 3'm where r > 1 and gcd (3, m) = 1. Then 
$(3n) = $(3"* 'm) = ф(3"*1)ф(т) = 3*3: (m) = 2-3" (m), 
and 
3ф(п) = 36 G*m) = 3:3" -$6(m) = 2-3'$(m) = p(n). 


We are also asked to show, conversely, that if $ (3n) = 3$ (n), then 3|n. This is 
the same as showing that if 3n then $ (3n) = 3ф (n). This is an immediate 
consequence of our solution to part (d). 


(d) If 3yn, then ged (3,n) = 1 so 
$n) = $(3) b(n) = 26(n) (#3¢(n)). 
Conversely, if ф(Зп) = 2(n), then ф(3п) = Зф(п) and so by part (с), 3n. 
(e) If n = 2, k = 1, then p(n) = $Q*) = 2*71 = 2 Conversely, suppose 


p(n) = 7 Then since $ (n) is an integer, we know that 2|n and we can write 


п = 2*N where k > 1 and N is odd. Then ф(п) = $(2*) $(N) = 2*- ! (N). 
But p(n) = 2 = 


5 2*-!N which gives Ф(№) = N. Now if N > 1 then 


1 
P2 


1 

во = ? | |: jep L) <x, and so we conclude that 
1 Pr 

N = 1 and n = 2*, as required. 


B page 141, Problems 7.2, number 13. 


‘Let n have prime factorization п = pp}? --- p*. By Theorem 6-1, any divisor 


d of n has the form 


d= pi'p?...pr, where 0 <a; € К, (for i= 1,2, ..., р). 


EXAMPLE 5 


SAQ 10 
SAQ 11 
SAQ 12 


SAQ 13 


SOLUTION 
TO SAQ 10 


SOLUTION 
TO SAQ 11 


SOLUTION 
TO SAQ 12 
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From the multiplicative nature of ф, 


oln) = фр) ф(р?)...ф(ру) 


апа 


pd) = ф(рү') 6(p7) ... (pr). 
To deduce that $(d)|ó(n), it suffices to show that $(p?)|ó(pF), where 
0 < a; < k; The case a; = 0 is easily dealt with, since then р“ = 1 and 
Q (p?) = 1. On the other hand, if a; > 0 then 
ФОР) = pr^! pe (p = 1) = p “(pf 


kia; 


(p — 1) = pi 
B page 141, Problems 7.2, number 10. 


If pi, pz, ..., p, are the distinct primes which divide m, then, since every prime 
dividing n also divides m, the same r primes are the ones which divide mn. 


Therefore 
1 1 
Жр, 
P2 Py 


ф(т) »[ E |! 
pi 


and 


(mn) m | 
Рі 
so ф(тп) = пф(т). 
Putting т = n we obtain p(n) = nọ (n). 
В page 140, Problems 7.2, number 1. 
B page 141, Problems 7.2, number 5. 


B page 142, Problems 7.2, number 17. 
[Note that the square brackets in line 2 of the question do not refer to the 
greatest integer function. | 


B page 142, Problems 7.2, number 19. 


1001 = 7:11:13, so 


(1001) = 1001(1 — 2) (1 — 45) (1 — 15) = 720. 
5040 = 24 -32-5 -7, so 
4(5040) = 5040(1 —3)(1 =- (1 23) 1 — 4) = 1152. 
36000 = 25-37-53, so 
(36000) = 36000(1 — 4)(1 — $) (1 —4) = 9600. B 
p(n) = b(2(2p — 1) = ф(2)ф(2р — 1) = (2p — 1) = 2p — 2, 


since 2p — 1 is prime. 
p(n + 2) = $4p) = $(4) b(p) = 2(p — 1) = 2p — 2, 
since p is prime. So ф(п) = ф(п + 2). E 
Let n = рер --- pF” so that 
b(n) = p p? pr (pi = 1) (ро — 1): (p, — 1). 


If ф(п) = 16, then (p; — 1)|16, so the only possible values for the p; are 2, 3, 5 
and 17. Furthermore, pF "|16 so either р; = 2, ог k; = 1. Therefore n must be of 
the form 2°3°5°174 where b, c and d are either 0 or 1. 


Suppose d = 1. Then n = 17m and $ (n) = $(17) $ (m) = 169 (m) so that 
o (m) = 1. Theorem 7-4 now implies that т = 1 or 2. It follows that 17 and 34 
are the only solutions involving the prime 17. 


SOLUTION 
TO SAQ 13 


5.4 


SAQ 14 
SAQ 15 
SAQ 16 
SAQ 17 


SOLUTION 
TO SAQ 14 
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We now assume d = 0, i.e. п = 243°5° where b and c are 0 or 1. Consider the 
case c = 1, so that 16 = p(n) = ф(2°3°) ф(5) = 49 2*3"), giving $(2^3^) = 4. For 
b = 0, this requires a = 3 and for b = 1, a = 2, giving respectively n = 40 and 

п = 60. Lastly, if c = 0, п = 273. For b = 1, 16 = p(n) = ф(2°) ф(3) = 2" giving 
а= 4 and n = 48. For b = 0, 16 = p(n) = $(2") = 2°! giving a = 5 and n = 32. 
So the only solutions of ф(п) = 16 are п = 17, 32, 34, 40, 48 and 60. 


Repeating the same arguments when ф(п) = 24, we see that the primes involved 
this time are 2, 3, 5, 7 and 13, so that n is of the form 2^3^5*7413*, Also, since 
5424, 7424 and 13/24 we must have that c, d and e are equal to 0 or 1. 
Looking each time at the largest prime involved, the case e — 1 leads to 
solutions 39, 52 and 78, the case d — 1 leads to solutions 35, 56, 70 and 84, 

с = 1 gives 45 and 90, and finally c = d = e = 0 gives solution 72. В 


o(p") = p" (p — 1). 
Therefore 


PP) = Ф(р“ 1) ф(р – 1), since gcd (p7 1,p — 1) = 1 
=p *(p — 1) b(p — 1). 
Now in Example 5 we showed that ф(п?) = пф(п), and so 
(p = 1) $(p — 1) = $((p — 1)?) 


giving 
PH) = p"?$(p—19) W 


Euler’s Theorem 


The remaining two sections of Chapter 7 of B investigate properties of the 
Q-function. We shall study only Section 7.3, which is concerned with Euler's 
generalization of Fermat's Little Theorem. 


COMMENT You may well recognize Euler's Theorem as being a particular 
case of a general result in group theory. Lagrange's Theorem (in group theory) 
tells us that if G is a finite group with order |G| and g is any element of G, then 
the order of g divides |G|. An immediate corollary of this is that g!^! = e, the 
identity element of G. Euler's Theorem is just this result for the case when G is 
the multiplicative group of units in Z,, whose order, as we remarked in the 
introduction to Section 5.3, is ф(п). To put the result in its right perspective, 
however, it must be remembered that Euler proved his result over a hundred 
years before the established birth of group theory. 


READ B Section 7.3, page 142 as far as page 145, line 10. 


B page 147, Problems 7.3, number 1. 

Find the final two digits of 13122. 

B page 147, Problems 7.3, number 3. 

B page 147, Problems 7.3, number 6. 

(a) It is sufficient to show that a?’ = a (mod p) for each of the primes p = 7, 13 
and 19, so in fact we need only use the Corollary to Fermat's Theorem, 


ie. a" = a (mod p). 
a?! = (a^a? = а?а? =a" = a (mod 7), 


аз? = (a? Pat! = aa" = a? = a (mod 13) 


and 


аз? = а!%а!% = aat! = ql? = a (mod 19). 
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(b) Similar to (a). We use the Corollary to Fermat’s Theorem to show that 
a? = a (mod p) for p = 2, 3, 5, 7 and 13. 


a? = (а?)ба = аба = (a° Pa = а?а = (a? = а? = a (mod 2), 


a!? = (а 


34 


а = ata = аа? = аа? = а? = a (mod 3), 


at? = (а?аз = а?аз = а? = a (mod 5), 


at? = а?а = аа = a? = a (mod 7), 


апа 
at? = a (mod 13). 
(c) In this case we do need Euler's Theorem, to deal with the factor 16. Since 


$ (16) = 8, we know that a? = 1 (mod 16) for any integer a such that 
gcd (a, 16) = 1, i.e. for any odd integer. So for any odd integer a we have 


a3 = (a®)*a = 1*a = a (mod 16). 


To deal with the remaining prime factors we revert to the Corollary to 
Fermat’s Theorem: 


а?? = (a3)! = a!! = (a? Ja? = aa? = аа? = à? = a (mod 3), 


a3 = (a°) a? = абаз = а?а = aa^ = а? = a(mod 5) 


and 


а3% = 211416 = qq16 = а!" = a (mod 17). 


This shows that a°? = a (mod 3۰5۰17) for all integers a, and 
а? = a (mod 3۰5۰16۰17) for any odd integer a. B 


SOLUTION We need to find the remainder on dividing 13!?? by 100. Since ф(100) = 40 and 
TO SAQ 15 gcd (13, 100) = 1, Euler’s Theorem shows that 134° = 1 (mod 100), and so 


13122 = (134°) - 13? = 13? = 69 (mod 100). Hence the final two digits of 1312 
аге 69. B 


SOLUTION Since gcd (т, п) = 1, we have n?" = 1 (mod m) and т = 1 (mod n). That is, 

TO SAQ 16 m|(n?™ — 1) and n|(m?™ — 1), giving mn|(n*” — 1) (т? — 1). In other words, 
mn|m? 9e — т) — п?) + 1, and since mn|m?n?™ this means that 
mn|( — m?*? — n?™ + 1). That is, m?? +n?” = 1 (поа тп). Wi 


SOLUTION (а) Substituting x = ba?” ! we have 
ИТ ах = aba®™~! = ba? = b (mod n). 
Thus x = ba? ~ is the (unique modulo n) solution of ax = b (mod n). 
(b) $(26) = 12, and so the solution of 3x = 5 (mod 26) is 
x = 5:311 = 5-(33)3-3? = 5-13-32 = 19 (mod 26). 
$ (40) = 16, and so the solution of 13x = 2 (mod 40) is 
x = 2-1315 =2-(137)7-13 = 2-97-13 =2-(9?)3-9-13 = 18:13 = 34 (mod 40). 
(49) = 42, and so the solution of 10x = 21 (mod 49) is 
x = 21-104! = 210: (102? = 210-22? = 14: (219) 
= 14-(—5)? = 7 (mod 49). Mi 


READ В page 145, line 11 to the end of Section 7.3. 


NOTES (i) Page 145, line 13 
The second proof of Euler's Theorem should be read for interest and as an 
exercise in manipulation. The previous, simpler, proof is preferred. 
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(ii) Page 146, line —3 
The Chinese Remainder Theorem is primarily concerned with the existence 
and uniqueness of a solution. However, here (as in the earlier proof on B 
page 87) the proof includes the construction of an actual simultaneous 
solution. But the solution produced will usually be inconveniently large 
and has to be reduced modulo n,n ...n,. For instance, consider the 
problem of solving simultaneously 


x =1(mod3), x = 3(mod 5), х =12(mod 14). 


In this case, №, = 70, №, = 42, N; = 15, ф(п,) = 2, (nz) = 4 and 
ф(пз) = 6. The solution given (by line —9, page 146) is 


x = 70? + 3-42* + 12-155 = 146027488 
which reduces modulo 210 to 208. 


So the method given here should be looked upon as theoretical; for 
practical purposes we revert to the technique of Unit 3. For the above 
problem: 


x = 12 (той 14) = x = 12,26, 40, 54, 68 (increasing by 
14 until 
division by 
5 leaves 
remainder 3) 

х=3 (mod5) = х = 68, 138,208 (increasing by 
5-14 = 70 
until division 
by 3 leaves 


x=1 (mod3) = x = 208. remainder 1) 


(iii) Page 147, line 9 
B remarks correctly that 7|111111. However, if you follow through the case 
n = 7, what is actually proved is that 7 divides 
10963 -1 109-4 
9 9 
The argument presents a general proof of the claimed result for any integer 


n relatively prime to 10. This result can be tightened when n is prime, as 
we see in SAQ 18. A 


ed (36 ones). 


SAQ 18 B page 147, Problems 7.3, number 7. 


SOLUTION Firstly, the result is certainly true for the prime 3 since by the divisibility test 3 
TO SAQ 18 divides every third term of the given sequence, namely, the terms consisting of 3 
ones, 6 ones, 9 ones etc. 


Suppose p is a prime other than 2, 3 or 5. Then by Euler's Theorem, since 
ged (10, p) = 1, 


109? = 1071 = 1 (mod p); 


that is, p divides 10^"! — 1 = 999...9, the number consisting of p — 1 nines. 
Now since р # 3, gcd (9, p) = 1 and by Euclid's Lemma, p divides a = 111... 1, 
the number consisting of p — 1 ones. 


Now the number consisting of 2(p — 1) ones is a + 107^! a, so it also is divisible 
by p. Similarly, we see that p divides every (p — 1)th term of the given sequence, 
namely the terms consisting of (p — 1) ones, 2(p — 1) ones, 3(p — 1) ones, and 
soon. Ё 


5:9 
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Summary 


In this unit we have studied some of the important functions of number theory, 
the central theme being the multiplicative functions. 


Section 5.1 was concerned with two functions, т and о. For any positive integer 
n, t(n) is the number of positive divisors of п, and o(n) is the sum of those 
divisors. Both functions were shown to be multiplicative and we obtained the 
following formulae in terms of the prime factorization n = pp? p: 


t(n) = (ky + UW (ky + 1): (k, + 1), 


Fatt _ 1 petii „+1 _ 1 
co(n) Pi .P2 „Pr 
р-1 p.-1 р. = 1 


Pursuing the o-function, we looked briefly at perfect numbers, which are positive 
integers n for which c (n) = 2n. We classified the even perfect numbers as being 
those integers of the form 2*^! (2* — 1), where 2* — 1 is prime. The primes of the 
form 2* — 1 are the Mersenne primes. 


In Section 5.2, we investigated properties of the greatest integer function and in 
particular used the greatest integer function to find the prime factorization of n! 


Section 5.3 introduced Euler's $-function, ф(п), which assigns to each positive 
integer n the number of elements in the set of least positive residues which are 
relatively prime to n, ф was shown to be multiplicative and we obtained the 


formula 
1 1 1 
1 1 {1 ——}, 
m 1 al A | _ 


where ру, p2,...,p, are the distinct primes dividing n. 


Section 5.4 presented Euler's generalization of Fermat's Little Theorem and 
looked at numerous applications. Euler’s Theorem states that if n is a positive 
integer and gcd (а, п) = 1, then а" = 1 (mod n). 


The main definitions and results encountered in this unit are listed in the 
corresponding entry in the Handbook. 
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UNIT 6 
THE QUADRATIC 
RECIPROCITY LAW 


6.0 


6.1 
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Introduction 


In Section 4.4 of Unit 4, pursuing applications of Fermat’s Little Theorem, we 
began a study of the solvability of polynomial congruences. We discovered 
similarities with the solutions of polynomial equations in that, provided we take 
congruences to a prime modulus, a polynomial congruence of degree n has at 
most n incongruent solutions. We now pick up this thread again, the central 
theme of this unit being the solvability of the general quadratic congruence. 


From the outset, we shall appreciate that the essential problem we must solve 
concerns the congruence 


x? = a (mod p). 


In other words, we need to know which elements in the arithmetic modulo p are 
the perfect squares. In the early units of the course, we solved such problems for 
several small values of p. For instance, we discovered that the square of any 
integer must take one of the forms 5k, 5k + 1 or 5k + 4. That is, in arithmetic 
modulo 5, the squares are 0, 1 and 4. It follows that the quadratic congruences 


x? =2(mod5) and х? =3(mod5) 
are each unsolvable. 


As the solution of quadratic equations presents little or no challenge, the 
thought of devoting a whole unit to quadratic congruences might well seem 
unexciting. Be patient! Investigation of when a quadratic congruence is solvable 
(without actually solving it) leads us through three of the outstanding results of 
number theory: Euler’s Criterion, Gauss’ Lemma and the Quadratic Reciprocity 
Law. The latter is a most fascinating result, described by Gauss as ‘the gem of 
higher arithmetic’. The consequences of these results are far-reaching in other 
areas of number theory. In this unit we shall use them to establish some more 
particular cases of Dirichlet’s Theorem (which we promised in Unit 2). There 
are numerous applications to the solution of Diophantine equations, some of 
which we shall meet in Unit 8. 


The unit is based on Chapter 9 of B. Section 6.1 looks at Euler’s Criterion, and 
introduces some new terminology. Section 6.2 presents the Legendre symbol (a 
most convenient shorthand for the new terminology) and investigates its 
properties. Section 6.3 proves Gauss’ Lemma and numerous corcllaries of it. 
Finally, Section 6.4 is concerned with the Quadratic Reciprocity Law and its 
applications, and is a tape-based section. 


In Chapter 9 of B the author gives several alternative proofs of results whilst 
pursuing a subsidiary line involving primitive roots. As we have chosen to omit 
primitive roots from our course, we have to steer a careful path through the 
chapter. 


Euler's Criterion 


In this section we obtain a first test, known as Euler's Criterion, for deciding 
whether or not the congruence х? = a (mod p) has a solution. In real terms, 
Euler's Criterion is not a very practical method to employ when the prime p is 
large, but nevertheless the test will prove crucial for subsequent theoretical 
considerations. 


We start, however, with a look at the problem of solving a general quadratic 
congruence, and introduce some new terminology. 


READ B Section 9.1, page 184 as far as page 186. line 19. 
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NOTES (i) Page 184, bottom 
The preceding argument has shown that the general quadratic congruence 
(1) is equivalent, after appropriate changes of variable, to congruence (2). 
The latter has solutions when the integer d is congruent to a ‘square’, 
modulo p. Therefore to solve the quadratic congruence we require a 
knowledge of which integers are ‘the squares, modulo p’. 


Notice that if yo is a solution of congruence (2), then to recover a solution 
of the equivalent congruence (1), via the equation y = 2ax + b, we have to 
solve the linear congruence 2ax = y, — b (mod p). But since ged (2a, p) = 1, 
this congruence has a unique solution (Corollary, B page 84). The 
numerical example in Note (ii) will illustrate the various stages of this 
argument. 


(ii) Page 185, line 14 
In this numerical example, B has taken a = 5, b 6, c = 2 and p = 13, 
and substituted these values into the various formulae on page 184. We 
shall look at the example in more detail. To solve the congruence 
5х2 — бх + 2 = 0 (mod 13), we first multiply through by 4а = 20 obtaining 


100x? — 120x + 40 = 0(mod 13). 


This allows us to ‘complete the square’: 
(10x — 6)? + 4 = 0 (mod 13). 
Adding 9 to each side, this becomes 
(10x — 6)? = 9 (mod 13). 
(At this stage in the accompanying theory we have d — 9 and we make the 
substitution y — 10x — 6.) 
This gives two possibilities, namely 
10x — 6 = 3 (mod 13) and 10x — 6 = —3 (mod 13). 
(Lagrange's Theorem shows that these two are the only possibilities.) These 
AR 10x = 9 (mod 13) and 10x = 3 (mod 13) 


and since gcd (10, 13) = 1, these linear congruences have unique solutions, 
which are respectively 


x =10(mod13) and x=12(mod13). A 


So the congruence ах? + bx + с =0(mod p) is solvable if and only if b? — 4ac is 
a quadratic residue of p. If this is the case, then the method given on B page 

184 (illustrated by the example in Note (ii)) will always lead systematically to 
the solutions. However, there are many other ways of solving quadratic 
congruences. For example, look again at the congruence 


5x? — 6x + 2 = 0 (mod 13). 


As a first step we can change the coefficients by multiplying through by any 
integer which is not a multiple of 13. Probably the simplest solution stems from 
multiplying by 8: 


40x? — 48x + 16 = х? + Ax +3 
= (х +2) — 1 =0(mod 13). 
Therefore.x + 2 = +1(mod 13), giving x = 10, 12 (mod 13), as before. 


We chose to multiply through by 8 simply because 8 is the multiplicative inverse 
modulo 13 of the coefficient 5 of x?, and therefore it converts the quadratic 
congruence to an equivalent one where the coefficient of x? is 1. The advantage 
of this is illustrated by the ease with which we completed the above solution. 


ко 
درا‎ 
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In our Solution to SAQ 1 we suggest various other ways of attacking quadratic 
congruences. There is also a point worth noting in the Solution to SAQ 2 about 
the number of quadratic residues of an odd prime p. 
SAQ 1 B page 189, Problems 9.1, number 1. 
, SAQ 2 B page 189, Problems 9.1, number 2. 
SOLUTION (а) х? + 7x + 10 = x? + 18x + 10(mod 11). (since 7 = 18 (mod 11)) 
HOSA = (х +9) —71 
= (x +9)” — 5 = 0 (mod 11). 


Therefore (х 4-9)? = 5 = 16 (mod 11), whence x +9 = +4 (mod 11), giving 
x = 6 (mod 11) or x = 9 (mod 11). 


Alternatively, the quadratic can be factorized: 
x? + 7x +10 = (х + 5)(х + 2) = 0 (mod 11). 


By Theorem 3-1 of B page 47, it follows that either x + 5 = 0 (mod 11) or 
X +2 = 0 (mod 11), giving x = 6,9 (mod 11), as before. 


(b) We first convert the coefficient of x? to 1. We can achieve this either by 
multiplying through by 9 or by cancelling common 3s: 


3x? + 9x +7 = 3x? + 9х + 33 (mod 13) (since 7 = 33 (mod 13)) 
= 3(x? + 3x + 11) = 0 (mod 13). 
Since 13/3, this is equivalent to 
x? +3x +11 = x? + 16x +11 
= (x +8)? — 53 
= (х +8)? — 1 = 0 (mod 13). 
Therefore x + 8 = + 1 (mod 13), giving x = 4,6 (mod 13). 
(c) Multiplying through by 14 (since 5-14 = 1 (mod 23)), 
70x? + 84x + 14 = x? — 8x + 14 
= (x — 4)? — 2 = 0 (mod 23). 


Therefore (x — 4? = 2 = 25 (mod 23), and so x - 42 +5 (mod 23) giving 
x = 9, 22 (mod 23). 


Alternatively, the quadratic again factorizes: 


5x? + 6x +1 = (5х + 1)(x + 1) = 0 (mod 23). 
So either 


X + 1 = 0 (mod 23) => x = 22 (mod 23) 
or 


5x + 1 = 0 (mod 23) => x = 9 (mod 23), 
giving solutions x = 9, 22 (mod 23). B 


SOLUTION (a) In Example 9-1, we determined the quadratic residues of 13 by evaluating 
TO SAQ 2 (modulo 13) the squares 17, 27,...,12?. In the same way, for апу odd prime 
p, the quadratic residues of p are given by the squares 12, 22,. .-,(p — 1). 
However, since r = (—r)? = (p — г)? (mod p), these p — 1 integers fall into 
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congruent pairs modulo p, namely 1? = (p — 1), 2? = (p — 2)?,. 


Бү 


Р—1\? 
ا‎ 
ER 


In fact, the integers listed above are incongruent modulo p, so that there are 
p-1 


2 
, and so a quadratic residue of p is congruent to one of 


precisely 


quadratic residues of p. To see this, suppose r? = s? (mod p), 


where 1 <r, s < P1 This means that p|r? — s?, and factorizing 

r? — s? = (r — s) (r + s), Euclid's Lemma gives p|(r + s) or pir = s). Now 
2<r+s<p—l1,and so pir + s is impossible. Moreover, |r — s| < p 
and so p|r — s can occur only if = s. 


(b) From part (a), the quadratic residues of 17 are congruent modulo 17 to the 
squares of the first eight natural numbers: 


12 = 1, 2? = 4, 32 = 9, 4 = 16, 5 = 8, 6 = 2, 72 = 15, 82 = 13. M 


Before reading about Euler’s Criterion, let us look at a couple of Examples. 


In SAQ 2 (b), we saw that 3 is a quadratic nonresidue of 17. Consider the linear 
congruence сх = 3 (mod 17), where c is any of the integers 1,2,...,16. Since 

gcd (с, 17) = 1, this congruence has a unique solution. The solution is certainly 
not x = 0 (mod 17), and neither is it x = c (mod 17), for then c? = 3 (mod 17) 
contradicts the fact that 3 is a quadratic nonresidue of 17. Hence, for each 
choice of c, the unique solution is one of the other integers in the set 

{1,2,..., 16}. It follows that these sixteen integers fall into eight pairs in such a 
way that the product of each pair is congruent to 3 (mod 17). To be precise, 


1:3 = 3 (mod 17), 14-16 = 3 (mod 17), 
2:10 = 3 (mod 17), 7:15 = 3 (mod 17), 
4:5 = 3 (mod 17), 12:13 = 3 (mod 17), 
6:9 = 3 (mod 17), 8-11 = 3 (mod 17). 
Multiplying together these eight congruences, we obtain 
16! = 38 (mod 17). 
But we know from Wilson's Theorem that 16! = — 1 (mod 17), and so 
38 = —1 (mod 17). 


The chosen integer 3 is not special in this example. If we select instead any 
other quadratic nonresidue of 17, say a, then we find that the sixteen integers 
from 1 to 16 inclusive fall into eight pairs whose product is congruent to 
a(mod 17). Multiplying together the resulting eight congruences gives 


a® = 16! = —1(mod 17). 


What is more, we can replace 17 by any odd prime p and discover the following 
result. 


If a is a quadratic nonresidue of p then a^^? = —1 (mod p). 


This is part of Euler's Criterion; the remaining part comes from a 
corresponding result when a is chosen to be a quadratic residue of p. 
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8 is a quadratic residue of 17. In SAQ 2 we found that the congruence 

x? = 8 (mod 17) has two solutions: x = 5 (mod 17) and x = 12 (mod 17). If we 
ignore 5 and 12 then, just as in Example 1, the remaining fourteen integers in the set 
11,2,..., 16} fall into seven pairs whose product is congruent to 8 modulo 17: 


1-8 = 8(mod 17), 9-16 = 8(mod 17), 
2۰4 = 8(mod 17), 10-11 = 8 (mod 17), 
3:14 = 8(mod 17), 13:15 = 8 (mod 17). 
6: 7 = 8(mod17), 


If we multiply together the left-hand sides of these congruences, the outcome is 
not 161, since the integers 5 and 12 are missing. To rectify this, and so that we 
can subsequently appeal to Wilson's Theorem, we include the additional 
congruence whose left-hand side is 5 times 12: 


5:12 =5-(—5) = —5? = —8 (mod 17). 


(It is no accident that the right-hand side turns out to be the negative of the 
quadratic residue under consideration.) 


Multiplying together all eight congruences, we have 
16! = 87: (—8) = — 88 (mod 17) 
and using Wilson's Theorem (16! = — 1 (mod 17)), we arrive at 
85 = 1 (mod 17). 


Again we can replace 17 by any odd prime p and 8 by any quadratic residue of 
р, say a, and the above argument generalizes to show 


a?~ 2 z1(modp) O 


The proof of Euler's Criterion in the next reading passage consists of the 
generalizations of our two examples. The statement of Euler’s Criterion in B is 
as follows: 


a is a quadratic residue of the odd prime p if and only if a'?~ 0/2 = 1 (mod p). 


We shall, however, prove the stronger formulation given as the Corollary on page 
187. 


READ В from the statement of the Corollary on page 187 to the end of Section 9.1. 


B page 189, Problems 9.1, number 8 part (a). 


When dealing with quadratic residues and nonresidues of p it is important to 
remember that we exclude from consideration the multiples of p. Here, since we 
are given that ab is a quadratic residue, it follows that pa and p/b, for 
otherwise p|ab, which is not the case. So neither a nor b is a multiple of p. 


Applying Euler's Criterion to the identity 
(aby»- 1/2 = ар 2 . pl 1/2 


we obtain the quadratic nature (residue or nonresidue) of ab in terms of that of 
a and b. For example, if a is a nonresidue and b a residue, then 


a 1/2 = ~], plr 2 = (ab)? 02 Ex (mod p) 


showing ab to be a nonresidue. Looking at the four possible cases, we draw up 
the following multiplication table: 
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x | residue nonresidue 
residue residue nonresidue 
nonresidue nonresidue residue 


which is determined essentially by the various products of 1 and —1. 


The table answers the question: residues arise only as the product of two 
integers of the same quadratic nature. 


Determine whether 5 is a quadratic residue or nonresidue for each of the primes 
p = 11, 13, and 17. 


If a is a quadratic residue of the odd prime p, prove that p — a is a quadratic 
residue of p if and only if p = 1 (mod4). 

B page 190, Problems 9.1, number 8 part (b). 

Lagrange's Theorem tells us that each of the polynomial congruences 

ХФ 0/2 — | = 0 (mod p) and x'"7 92 + 1 = 0 (mod p) has at most 2 3 > 


incongruent solutions. However, 
(х D? — 1) (x79? 4. 1) = xP“! — 1 =0(modp) 


has exactly p — 1 solutions, since x^! = 1 (mod p) for x = 1,2,...,p — 1, by 
Fermat's Theorem. 


By showing that each quadratic residue of p is a solution of 
x” 0? —j = 0 (mod p), use the above information to give an alternative proof 
of Euler's Criterion. 


501-0? = 5° = (25)? -5 = 32.5 = 45 = 1 (mod 11), 
so 5 is a quadratic residue of 11. 
503717? — 56 = 253 = (—1)? = —1 (mod 13), 
so 5 is a quadratic nonresidue of 13. 
507-0 = 58 = 254 = 84 = 64? = 132 = 169 = – 1 (mod 17), 


so 5 is a quadratic nonresidue of 17. Bill 


p — а= —a(mod p). Therefore 


(р а)? 02 = ( ayr- 0/2 = (-1) 2-0? а 02 


= (—1)""? (mod p), 
since a?~"? = 1 (mod p) by Euler's Criterion, a being a quadratic residue. 
If p = 1 (mod 4), say p = 4k + 1, then 
(- 197? = (19% 21, 
in which case p — a is a quadratic residue of p. 
On the other hand, if p = 3 (mod 4), say р = 4k + 3, then 
( 1)7 1/2 ( jja 1, 


so that p — a is a quadratic nonresidue of p. W 


Consider the congruence ax? = b (mod p), where a and b are of the same 
quadratic nature modulo p. Being either a quadratic residue or nonresidue, 
a £ 0 (mod p) and so (as we have seen several times before, most recently in the 
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proof of Euler’s Criterion) a has a multiplicative inverse a’, that is, an integer a’ 
such that aa’ = 1 (mod p). Now 1 = 1? is certainly a residue of every prime. As 
the product aa’ is a residue of p it follows from Example 3 that a’ and a (and 
therefore b also) are of the same quadratic nature. 


Multiplying through the given congruence by a’ (obtaining an equivalent 
congruence since gcd (a’, p) = 1), we have 


аах? = x? = a'b (mod p). 


But as a’ and b are of the same quadratic nature, their product is known to be 
a residue, and so х? = a'b (mod p) certainly has a solution. Bl 


SOLUTION From the information given, each of the polynomial congruences 

TO SAQ 6 x?" "?—1-0 (mod p) and x??? +1 = 0 (mod p) has exactly 4(p —1) 
solutions. We have seen in the Solution to SAQ 2 that there are exactly 
3(p — 1) quadratic residues of p. If we can show that each of the 4(p — 1) 
quadratic residues is a solution of х- 1/2 —1 = 0 (mod p), it will follow 
that the remaining }(p — 1) non-zero least positive residues, which are 
the quadratic nonresidues, must each be a solution of х7? + 1 = 0 (mod p). 


Let x be a quadratic residue of p. Then x = a? (mod p), for some integer a with 
gcd (a, p) = 1. But then 


x97? = (42) 02 = ар-1 = 1 (mod p), 
by Fermat's Theorem, and so х 1? — | = 0 (mod p). 
It follows that a is a quadratic residue or nonresidue of p if and only if, respectively, 
а? 02 — 1 =0(тойр) ог a*-9? + 1 = 0(modp). 


This is Euler's Criterion. Bi 


6.2 The Legendre Symbol 


Throughout the previous section, you will have noticed the large number of times 
we wrote down statements such as ‘a is a quadratic (non)residue of p. You 
probably also noticed that from time to time we dropped the adjective 
‘quadratic’ and referred simply to ‘residues’ and ‘nonresidues’. It would be most 
convenient to have a shorthand notation for statements such as the above. To 
this end we shall introduce a very useful symbolism first used by Legendre. 


The section pursues our investigation of properties of quadratic residues using 
the new notation. 


READB from the beginning of Section 9.2, page 190, as far as the statement of | 
Theorem 9-4 on page 194. 


NOTES (i) Page 191, line —4 
The earlier example referred to is Example 9-1 on page 186. 
(ii) Page 192, line 16 
The proof given of Theorem 9-2 (1) is somewhat verbose. You might find 
the following argument clearer. 


(a/p) = 1<== а = х? (mod p), for some integer x 
«— b = х? (mod p), since a = b (mod p) 
<=> (b/p) = 1. 
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(ili) Page 193, line 17 
In Example 9-4, the alert reader might spot that —38 = 1 (mod 13), so that 
the congruence x? = —38 (mod 13) has solutions x = + 1 (mod 13). 
However, the example is intended to demonstrate parts of Theorem 9-2 in 
action and should be read in this spirit. The result itself is unimportant, 
but make sure you follow the steps taken in the argument. 


(iv) Page 194, line 14 
We proved Theorem 9-4 in our Solution to SAQ 2. There are 4(p — 1) 
quadratic residues of p, and they are congruent modulo p to the squares 


22... 


. The remaining (р — 1) incongruent numbers which 


are relatively prime to p are the quadratic nonresidues. 


p-1 
In the sum У (a/p), we are counting +1 for each of the 1(р — 1) residues 
=1 


and —1 for each of the ¿(p — 1) nonresidues. Clearly the total is zero. А 


EXAMPLE 4 As an example of the use of the properties of Theorem 9-2, let us find all the 


SAQ 7 


SAQ 8 


SAQ 9 


quadratic residues of 37. We know that there are eighteen of them, and we 
could, of course, find them by evaluating 17, 2,..., 18 modulo 37, but we hope 
we can reduce the amount of computation. 


For a start, we know that 1, 4, 9, 16, 25 and 36 are residues of 37 (this by (2)). 
As 36 = —1(mod 37), (1) gives (— 1/37) = 1, a fact that we could have 
established alternatively from (5). Using this together with (4) and (1), 


1 = (— 1/37) (4/37) = (— 4/37) = (33/37), 


showing that 33 is also a residue. In the same way, —9 = 28, —16 = 21 and 
—25 = 12 are residues. Having found that 12 and 28 are residues we can bring 
property (6) into play: 


1 = (12/37) = (2? - 3/37) = (3/37) 
and 

1 = (28/37) = (22 - 7/37) = (7/37). 
So 3 and 7 are residues, as too are —3 = 34 and —7 = 30. 


We have now found a total of 14 so there are four more outstanding. Using 
properties (4) and (1): 


1 = (4/37) (12/37) = (48/37) = (11/37) 

giving 11 and —11 = 26 as two further residues. 
Finally, using property (6) in a different way: 

(3/37) = 1 == (3? - 3/37) = (27/37) = 1 
giving 27 and —27 = 10 as the final two residues. 
The complete list of residues of 37 is 

1, 3, 4, 7, 9, 10, 11, 12, 16, 21, 25, 26, 27, 28, 30, 33, 34, 36. 

Using only the properties of Theorem 9-2, evaluate the following Legendre 
symbols. 
(a) (3/11), (Ы) (19/17, (е) (40/11). 


Show that a is a quadratic residue of 31 if and only if 2a is also a quadratic 
residue. Use this fact to obtain all quadratic residues of 31. 


B page 203, Problems 9.2, number 12. 
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(i) (3/11) = (3-22/11) (by (6)) 
= (12/11) 
= (1/11) (by (1)) 
=1 (by (5)) 
(ii) (19/17) = (36/17) (by (1)) 
=1 (by (2)) 
(ш) (40/11) = (4/11) (10/11) (Бу (4)) 
= (10/11) (by Q)) 
= (—1/11) (by (1)) 
= (–1)5 (by (5)) 
=-1 Ш 
(a/31) = (87a/31) (Бу (6)) 
= (64a/31) 


= (2a/31) (by (1)) 


So if we double (modulo 31) any residue of 31 the outcome is another residue. 
Starting from residue 1, doubling gives further residues 2, 4, 8 and 16. Starting 
from 9 = 37, doubling gives 18, 36 = 5, 10 and 20. Starting from 25 = 52, 
doubling gives 50 = 19, 38 = 7, 14 and 28. We have now found fifteen residues 
the correct number for the prime 31. They are 


1, 2, 4, 5, 7, 8, 9, 10, 14, 16, 18, 19, 20, 25, 28. B 
(p-1)2 
In the sum У, (a/p), we are examining the set of integers {1,2,...,4(p — 1)} 
a=1 
and counting | for each residue and —1 for each nonresidue. The total will be 
zero provided there are the same number of residues as nonresidues in this set. 


» 


In SAQ 4, we showed that if p = 1 (mod 4) then (a/p) = (p — a/p). Hence if a is 
a residue in (1,2,...,3(p — 1)} then p — a is a residue in {4(p + 1),...,p — 1}, 
from which it follows that exactly one half of the residues are to be found in 
{1,2,...,2(p — 1)}. This set therefore contains }(p — 1) residues and 1(p — 1) 
nonresidues. B 


Gauss" Lemma 


In this section we shall prove Gauss’ Lemma and look at some consequences. 
Gauss' Lemma gives another test for whether or not an integer is a quadratic 
residue of a prime, and, like Euler's Criterion, whilst not being a very 
practicable test when large numbers are involved, it is an important theoretical 
result. 


Before embarking on a proof of Gauss’ Lemma we shall look at an example 
which should clarify the statement of the result and, more important, will 
illustrate the line to be adopted in the general proof. 


READ B the statement of Gauss Lemma, Theorem 9-5 on page 195. 


Let p = 19 and a = 11. The set S in Gauss’ Lemma consists of the first 
(р — 1) = 9 multiples of 11: 


S = (11, 22, 33, 44, 55, 66, 77, 88, 99}. 
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Replacing each member of S by its remainder on division by 19, that is, by the 
least positive residue of 19 to which it is congruent, we get 


S' = (11, 3, 14, 6, 17, 9, 1, 12, 4}. 


There аге five members of S’ which are less than 3р and four members, namely 
11, 14, 17 and 12, which exceed 3p. Hence, according to Gauss’ Lemma, 


(11/19) = (-15 =1 


and 11 is a quadratic residue of 19. 


Suppose now that we replace the four members of S’ which exceed 4p by their 
negatives, achieved up to congruence by subtraction from 19. This results in the 
new set 


S” = {8, 3, 5, 6, 2, 9, 1, 7, 4}. 


Notice that 5” consists of the first nine natural numbers. This is no accident! 
The crux of the proof of Gauss’ Lemma lies in showing that, for any prime p 
and relatively prime integer a, the set S", constructed from multiples of a in this 
way, consists of the first 5(р — 1) natural numbers. It is a short step from here 
to completion of the proof. 


Returning to our specific example, notice that the product of the elements in S" 
is 9!, whilst the product of the elements in S is 11°-9! Now corresponding 
elements of the sets S and $” are congruent modulo 19, with the exception of 
the four sign changes involved in getting from S$’ to S”. It follows that 

2290 =f TW 
giving 117۰9! = (—1)*9! (mod 19) 

11° = 1 (mod 19), since ged (19,9!) = 1. 

Euler’s Criterion now confirms that 11 is a quadratic residue of 19. 


The proof of Gauss’ Lemma retraces this example in generality, for arbitrary 
odd prime p and any integer a with gcd (a, p) = 1. 


READ В page 195, line —7 to page 198, line —7. 


(i) Page 195, bottom 
The set of residues (ri,r2,..., м, S1,...,S,} corresponds to the set 5' of 
Example 5 (although the members are presented in a different order). 
ri Ро, sys D — S15+++»P — Sn} corresponds to our set S”. 


(ii) Page 196, line 15 
The use of square brackets here, and two further instances in the next eight 
lines, is accidental. They do not denote the greatest integer function. 


(11) Page 198, top 
We have found that the number of integers in the set {2,4,6,...,p — 1} 
which exceed 2р is ¿(p — 1) — [4p]. To apply Gauss’ Lemma we need to 
raise —1 to this power, and so what we really need to know now is whether 
(р — 1) — [р] is even or odd. To achieve this we have to look at the four 
forms for an odd prime modulo 8. It is not good enough, for instance, to 
work with the two forms 4k + 1 and 4k + 3, for if р = 4k + 1 then 


1р — 1) – Bp] = 2k — [k +4] = k 


and this can still be either even or odd. A 


SAQ 10 В page 201, Problems 9.2, number 1 parts (а), (c) and (e). 


SAQ 11 


B page 202, Problems 9.2, numbers 4 and 5. 


SAQ 12 В page 209, Problems 9.3, number 4. 
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As in the proof of Gauss’ Lemma, to determine (a/p) we form the set S 
consisting of the first 3(p — 1) positive multiples of a, reduce each to its least 
positive residue modulo p, and count how many resulting members exceed p/2. 
This gives the number n in Gauss’ Lemma. 


(а) a=8,p=11. 
S = (8, 16, 24, 32, 40} = (8, 5, 2, 10, 7). 
n= 3, so (8/11) = )-1( = –1. 
(c) a=5,p = 19. 
S = {5, 10, 15, 20, 25, 30, 35, 40, 45} | 
= {5, 10, 15, 1, 6, 11, 16, 2, 7}. 
n=4, so (5/19) =(—1)* = +1. 
(e) a= 6, p = 31. 
S = {6, 12, 18, 24, 30, 36, 42, 48, 54, 60, 66, 72, 78, 84, 90} 
= {6, 12, 18, 24, 30, 5, 11, 17, 23, 29, 4, 10, 16, 22, 28}. 
n=9, so (6/31) = )-1 = —1. B 


4. By Theorem 9-6, (2/p) = 1 for p = 7 (mod 8). By Theorem 9-2 (3), 
(2/p) = 2”? (mod p). Combining these results gives, for р =7 (mod 8), 
20-102 21 (mod p), or p|2- 9/2 — 1, 

5. The numbers in question, M, — 2" — 1, are the Mersenne numbers which we 
mentioned in Unit 5. If n is composite, M, is known to be composite, but for 
p prime M, might or might not be prime. The above result uncovers many 
composite values of M,, for it tells us that if p = 7 (mod 8) is a prime such 
that 5 (p — 1) is also prime, then М (p-1y2 İS composite (being divisible by p). 


For p = 23, (р — 1) = 11 is also prime and so 23|M ||. 
Similarly, M23, Маз, M131, M179, M183, M239 and M55, are shown to be 
composite by applying the above result to the respective primes 47, 167, 263, 
359, 367, 479 and 503. BH 
To apply Gauss’ Lemma we look at the set 
S={-2, -4, -6..., -(p- 1). 


Adding p to each member of S (and reversing the order), we obtain the 
congruent set, which are positive least residues: 


S — (1,35... pp 2) 


We need to count how many of these integers exceed p/2. To do this we look 
separately at the four possible forms of p modulo 8. 


G) Ifp = 8k +1, then the members of S' which exceed p/2 are 
4k + 1, 4k + 3,...,8k — 1, 
of which there are 2k. In this case, (—2/p) = (—1)?* = 1. 
(ii) If p = 8k + 3, then the corresponding list is 
4k + 3, 4k +5,...,8k +1, 
of which there are again 2k, giving (—2/p) = 1. 


(ii) If p = 8k +5, the list is 
4k + 3, 4k + 5,...,8k + 3, 


of which there are 2k + 1. In this case (—2/p) = (—1)**! = —1. 
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(iv) If p = 8k +7, the list is 
4k + 5, 4k +7,...,8k +5, 
of which there are again 2k + 1, giving (—2/p) = —1. 
Alternatively 
(—2/p) = (—1/p): (2/p) (Бу Theorem 9-2 (4)) 
= (– 1) 0/2. ) 1(2 0/8 (by Theorems 9-2 (5) and 9-6 (Corollary)) 


= (— ту 0/2)+ (9? = 1)/8) 


= (- 1)? *4»-5)8. 


Therefore (—2/p) = 1, provided р? + 4p — 5 = 0 (mod 16), which is the case 
when р = 8k + 1 or p = 8k + 3. On the other hand, (—2/p) = —1 provided 
p? + 4p — 5 = 8 (mod 16), which is the case when p = 8k + 5 огр = 8k + 7. W 


READ B page 199, line —6 to page 200, line 8. 


Page 200, line 5 
N = (pipi: p — 2 = 2Q psp; p, — 1). 


As the expression in curly brackets is odd, N is equal to twice a product of odd 
primes. If each of these odd prime divisors of N were congruent to 1 modulo 8, 
then so too would be their product, forcing N = 2(8a + 1) = 16a 4- 2. This 
cannot be the case since N is of the form 16а — 2. А 


For any given integer n > 2, let N =2 (n!) — 1, and let p be a prime divisor of 
N. By actually finding a solution of the congruence x? = 2 (mod p), show that 
(2/p) = 1. Hence, using Theorem 9-6, deduce that N must have a prime divisor 
of the form 8k + 7. 


Why does this argument prove that there are infinitely many primes of the form 
8k + 7? 


If p is a prime divisor of N = 2(n!)? — 1, then p is odd and 2(n!) = 1 (mod p). 
Multiplying by 2 gives (2(n!))* = 2 (mod p), showing that (2/p) = 1. According 
to Theorem 9-6, either p = 1 (mod 8) or p = 7 (mod 8). 


Now, if every prime divisor of N is of the form 8k + 1, then so is N. But N is of 
the form 8k — 1 (since 8|2(n!)’, for n > 2). Therefore N must have at least one 
prime divisor of the form 8k + 7. 


Suppose there are only finitely many primes of the form 8k + 7, and let n be the 
largest of them. Then 2(n!? — 1 is divisible by a prime p = 8k + 7 and, as 

p < n, pln! Therefore p|2(n!)?. This is a clear contradiction since it forces 

pi Ш 


Quadratic Reciprocity 


We now turn to what is undisputably one of the most fascinating results in the 
whole of mathematics—Gauss’ Quadratic Reciprocity Law. Amongst its 
numerous applications, the law provides us with yet another test for when an 
integer a is a quadratic residue of a prime, but, unlike its predecessors, this test 
is simple to employ even when the integers in question are very large. 


As you have already read, there are many different proofs of the Quadratic 
Reciprocity Law. We have elected to present here our favourite proof, which is 
slightly different from the proof given in B. If you wish to study the proof in B, 
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either instead of or in addition to our proof, then you will need to read first the 
Lemma on B pages 200 and 201. 


READB Section 9.3, page 203, as far as the statement of the 
Quadratic Reciprocity Law on page 204. 
START THE TAPE WHEN YOU ARE READY 
][-— 7-3 


FRAME 1 Gauss’ Quadratic Reciprocity Law 


If p and q are distinct odd primes then 
(p/a): (a/p) = (—1)e- 99-09%, 
that is 
(р/а) = (q/p) if either p or q is congruent to 1 modulo 4 
(р/а) = —(a/p) if p = q = 3 (mod 4). 


FRAME 2 Example 6 
Is 31 a quadratic residue of 73? 
(31/73) = (73/31), since 73 = 1 (mod 4) 
= (11/31), since 73 = 11 (mod 31) 
= —(31/11), since 11 = 31 = 3 (mod 4) 
—(9/11), since 31 = 9 (mod 11) 
= —1, since (9/11) = (32/11) = 1. 


FRAMES Restating the problem via Gauss’ Lemma 
Let m be the number of integers in the set 
1.24, ..., (p — 1)q} 
which have least positive residue modulo p exceeding p/2. 
Let n be the number of integers in the set 
{р,2р,...,3(9 — 1)p} 
which have least positive residue modulo q exceeding q/2. 
Gauss Lemma: (q/p) = (—1)"; (p/q) = (—1)". 
Now (p/q) = (4/p) => (p/q): (a/p) = (—1)"*" = 1. 


So to prove the Quadratic Reciprocity Law: 


Show that m + n is even except when 


P = q = 3 (mod 4). 
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FRAME 4 Diagram 


FRAME 5 


FRAME 6 


Y(o.$), 2 C(3) 
AB 
Е (0, 
|р 
° 40,0) x (5. о) 


OXCY is the rectangle with vertices (0,0), (p/2,0), (p/2, q/2) and (0, q/2), 
respectively. 


OC is a diagonal and AB and ED are parallel to OC passing respectively 
through (3, 0) and (0,4). 


Equations: 
А =a. _ 1 е „ыа 
ОС: Py AB: y= к= 2), ED: уке ete 


Plan of the proof 

Lattice points in OXCY are points with coordinates (a,b) where a and b are 
both positive integers. 

We shall show: 


(1) The number of lattice points inside the hexagon OABCDE is even, except 
when p = = 3 (mod 4). 

(2) The number of lattice points in quadrilateral OABC is n. 

(3) The number of lattice points in quadrilateral OCDE is m. 

(4) None of the lattice points inside OABCDE lie on OC. 


Outline proof of (1) 

The lattice points in OABCDE fall naturally into pairs, because if (a,b) is such a 

1 ud 1 
"29 


point, then so too is р — 2 These points are on the straight 


Р ag , 1 1 
line through, and equidistant from, the point E roe. fr 


41). (Ѕее Егате 4.) 


1 
ана аге ће 
1 11. i 
. Hence, if r, ч ы is a lattice 


There is one exception to this pairing. (a,b) and 


; 1 +1 
same point when a р м and b = 4 


point, which will be the case precisely when p = q = 3 (mod 4), this exceptional 
point is paired with itself. In this case alone the total number of lattice points in 
OABCDE will be odd. 
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Detailed proof of (1) 


The interior of the hexagon OABCDE consists of 


points (x, у) satisfying the 
constraints 


4 


р 
# 0 s 
А257. MATO 


q (See Frame 4.) 


1 4 
Ум Psy 
p 


q 
у> 2х -> 
р 


2° 2p" 
It is our contention that if x = 4, y = b satisfies these inequalities then so too 
+1 4+1 


гаје 5—3, 


does х =? 


, 1 
(since ps 


— a is an integer). 
Similarly, 


1 q«i 
i RAN uM =bx 


qti 


4+1 
2 


q 


Y 


q 
Os Baas 


=0 < b< 


q q 
Ь>-а——- 
р 2р 


апа 


Proof of (2) 


Objective: To show there is a one-one correspondence between the lattice 


; 3 — 1 ; 
points in OABC and the numbers in the set (p, 2p, 3p... 57 — p) which have 
least positive residue modulo q exceeding 5. 


If (а, Б) is a lattice point in OABC then 
(i) 
and 


(i) 
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Multiplying (i) by p and adding q — qa gives 
1 <р qa- 1)< q. (iii) 


Now bp = bp — q(a — 1) (mod q), and so (iii) says that bp has least positive 


residue modulo q exceeding 2 


— 1 
The argument reverses. Choosing b from the set 11.2 „ә (so that (ii) 


holds), if bp has least positive residue modulo q exceeding 3 then there exists an 


integer a such that (iii), and hence (i) holds. (a, b) is therefore a lattice point in 
OABC. В 


Proof of (3) 


The proof of (3) is almost identical to the proof of (2), with the roles of p and q 
interchanged. 


SAQ 14 
Prove (3). 


Solution to SAQ 14 
If (a, b) is a lattice point in OCDE then 


1 J 
Traba ағы (i) 
р р 2 
and 
р js 
4 i 
0<а<5 (1) 


Multiplying (i) by —p gives 
—qa> —bp > —qa -£ 
and adding qa + p gives 
p > qa — (b — Dp > 5. (iii) 


(iii) shows that qa has least positive residue modulo p exceeding p/2. 
TU — 1 
Conversely, if a is chosen from the set L 2, MM 251) so that ад has least 


positive residue modulo p exceeding 5 then there is an integer b such that (iii), 


and equivalently (i), holds. (a, b) is then a lattice pointin OCDE. W 


Proof of (4) 
If (a, b) is a lattice point on OC then b — a. Hence p|qa, and since 


gcd (p,q) = 1, pla. Inside OABCDE we have 0 < a < 5, and so pla is 


impossible. Mi 
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Example 7 
Does the quadratic congruence 


3x? — 8x — 5 = 0 (mod 139) 
have any solutions? 
b? — 4ас = (C8? — 4-3-( 5) = 124 


The congruence has solutions provided that 124 is a quadratic residue of 139. 


(124/139) = (4/139) - (31/139) 


= 1- (31/139) (since 4 = 22) 

= —(139/31) (since 139 = 31 = 3 (mod 4)) 
= —(15/31) (since 139 = 15 (mod 31)) 

= — (3/31) (5/31) 

= (31/3) (5/31) (since 31 = 3 (mod 4)) 

= (31/3) (31/5) (since 5 = 1 (mod 4)) 

= (1/3) (1/5) 


=1. 


The congruence has solutions. 


For more examples, 


end of Example 9-6 on page 208. 

END OF TAPE 
(a) Evaluate the following Legendre symbols: 

(i) (—219/383), (ii) (461/773). 

(383, 461 and 773 are primes.) 
(b) Are the following congruences solvable? 

(1) 3x? + 6x + 5 = 0 (mod 89) 

(ii) 2x? + 5x — 9 = 0 (mod 1111). 


If p and q are distinct odd primes, prove that (—p/q) = (—q/p) if and only if 
p =q = 1 (mod 4). 


B page 209, Problems 9.3, number 5. 
B page 210, Problems 9.3, number 8. 


(a) Explanations of the steps refer to properties (1), (2), (4) and (5) of the 
Legendre symbol (Theorem 9-2), Theorem 9-9 (GQRL) and Theorem 9-6. 


G) (—219/383) = (164/383) (by (1)) 
= (4/383) - (41/383) (by (4)) 
= 1- (383/41) (by (2) and GORL) 


= (14/41) = (2/41): (7/41) (by (1) and (4)) 
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= 1: (7/41) (by 9-6) 
= (41/7) (by GQRL) 
= (6/7) = (2/7) - (3/7) (by (1) and (4)) 
= 1-6/7) (by (9-6) 
= —(7/3) (by GQRL: note 7 = 3(mod4)) 
= —(1/3)= -1. (by (1) and (5)) 
(ii) (461/773) = (773/461) = (312/461) (by GQRL and (1)) 
= (4/461) - (2/461) - (3/461) - (13/461) (12 = 2*-3-13) 
= 1-—1-(461/3)- (461/13) (by (2), 9-6 and 
GQRL) 
= —1-(2/3)- (6/13) (by (1)) 
= —1: —1 - (2/13): (3/13) (by 9-6 and (4)) 
= —1-(3/13) 
(13/3) (1/3) 1. (by GQRL and (1)) 


(i) 


and 


3x? + 6x + 5 = 0 (mod 89) is solvable provided 6? — 4۰3۰5 = —24 isa 
quadratic residue of 89. 


(— 24/89) = (65/89) = (5/89) - (13/89) 
= (89/5) - (89/13) 
= (4/5): (11/13) 
= (13/11) = (2/11) 2 —1 = No solutions. 
1111 = 11-101, and so the congruence 2x? + 5x — 9 = 0 (mod 1111) is 
solvable if and only if both the congruences 2x? + 5x — 9 = 0 (mod 11) 


and 2x? + 5x — 9 = 0 (mod 101) are solvable. This will be the case if 
5? — 4:2-(—9) = 97 is a quadratic residue of both 11 and 101. 


(97/11) = (9/11) = (32/11) 2 1 
(97/101) = (101/97) = (4/97) = (22/97) =1 


So the congruence has solutions. W 


p/a) = (— 1/4) (р/а) = (— 1)" P (p/q) (by Theorem 9-2 (5)) 


(—4/р) = (—1/p) (a/p) = (—1)** (дур). 


Therefore (— p/q) = (— q/p) if either 
@) (—1)«7? = (— 1)? and (p/q) = (q/p), or 
(i) (1) = —(—1)+Ф—1) and (р/а) = —(а/р). 


Taking case (i), (— 1)" = (—1)#”~” when p = q (mod 4), whilst (p/q) = (q/p) 
except when р = q = 3 (mod 4). We are left with the single solution 


р = = 1 (mod 4). 


Turning to (ii), (p/q) = —(q/p) only when p = = 3 (mod 4), but this possibility 
contradicts (—1)54 ^9 = —(—1)i(- 9, 


We conclude that (—p/q) = (—q/p) if and only if p = q = 1 (той4). Ш 
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SOLUTION (а) (—3/p) = (—1/p): G/p) 
TO SAQ 17 1-(p/3) 
Cd ee 
= (p/3) 


m ifp = 1 (mod3) 
1-1 ifp =2(mod 3), 


ifp = 1 (mod 4) 
ifp = 3 (mod 4) 


(by Theorem 9-2 (5) and 
GQRL) 


COMMENT This is the same as the given answer, for since a number of the 
form 6k + 4 is not prime, the set of primes which are congruent to 1 modulo 3 is 
precisely the set of primes which are congruent to 1 modulo 6. Similarly, if p is a 
prime, then p = 2 (mod 3) if and only if p = 5 (mod 6). The result of (b) also 
shows that there are infinitely many primes of the form 3k + 1. 


(b) Suppose pı, p2,..., р, are the only primes of form 6k + 1 and consider 


N = (2p1p2 py + 3. 


As N is odd it has an odd prime divisor, say p. Then 

Qpip; ... p, = —3 (mod p), from which (—3/p) = 1 and so p = 1 (mod 6) 
by part (a). But this means that p must be one of P1:P2:۰.., P, Whence 

PIN and p|(2p,p, ... р,)?, giving p|3. This implies that p = 3, contradicting 
p = 1 (mod 6). Hence there are infinitely many primes of the form 


6k+1. E 


SOLUTION (a) (5/p) = (p/5). The quadratic residues of 5 are 1 and 4, so 


TO SAQ 18 


COMMENT 
(b) (6/p) = (2/p)- 


(6/р) = 1 = 


(5/ (=1 = م‎ =1 or 4 (mod 5). 


Again this is a more compact form of the answer given in B. 
(3/p), so 


(2/р) = (3/р) = 1 
OR 
(2/p) = (3/р) = —1 
p=+1(mod8) AND p= +1 (mod 12) 
OR 
p= +3(mod8) AND p= +5(mod 12) 


p = 1 or 23 (mod 24) 
OR 
р = 5 or 19(mod 24) 


<= p = 1, 5, 19 or 23 (mod 24). 


(©) (ip) = iv 7) 


ifp = 1 (mod 4) 


—(p/T) ifp = 3 (mod 4) 


The quadratic residues of 7 are 1, 2 and 4, hence 


(7/р) = 1— 


چس 
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p=1(mod4) AND p=1,2 or 4(mod7) 
OR 
p=3(mod4) AND p33, 5 or 6(mod7) 
р = 1, 9 or 25 (mod 28) 
OR 
р = 3, 19 or 27 (mod 28) 
P = 1, 3, 9, 19, 25 or 27 (mod 28). M 
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6.5 Summary 


The central theme of this unit concerns the problem of determining whether or 
not a quadratic congruence with prime modulus has solutions. The answer to 
this problem hinges around our ability to decide whether a given integer is a 
quadratic residue (i.e. congruent to a square) or nonresidue of the prime 
modulus. 


In Section 6.1 we began our investigation of quadratic residues and showed that 
an odd prime p has $(p — 1) residues and 4(p — 1) nonresidues. We obtained a 
first general test for residues, Euler’s Criterion, which says that a is a quadratic 
residue of p if and only if 


a 92 = 1 (mod p). 


Section 6.2 pursued the development of the theory of quadratic residues, 
incorporating a new notation—the Legendre symbol: 


яў 1 ifais a quadratic residue of p, 
Wo] 4 ais a quadratic nonresidue of p. 


Theorem 9-2 (B page 192) contains the important properties of the Legendre 
symbol. 


In Section 6.3 we gave a proof of Gauss’ Lemma (B page 195), which is 
another test for the quadratic nature of an integer modulo p. Amongst several 
applications of Gauss’ Lemma we proved some more special cases of Dirichlet’s 
Theorem, that certain arithmetic progressions contain infinitely many primes. 


The climax of the unit was Gauss’ Quadratic Reciprocity Law (B page 204), 
which we proved in Section 6.4. This result, coupled with properties of the 
Legendre symbol, provides us with a relatively simple technique for deciding the 
solvability of any quadratic congruence with prime modulus. 


The main definitions and results encountered in this unit are listed in the 
corresponding entry in the Handbook. 
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Introduction 


The most obvious way of solving the quadratic equation x? — x — 1 = 0 is to 
use the well-known formula. Suppose, however, that we decide to solve the 
equation using an iterative method. The equation can be rearranged as 


1 " А Я r 
x = 1 + — so we could find successive approximations, x,, to a root by using the 
56 


formula 1 
Mati = 1+ —. 
Xn 


If we take x, = 1 as our initial approximation, we obtain the following: 


Ad m1 
1 I1 2 
= 1S 462-2 
s r E ERE 
1 3 
xum 1+1 5 
"а 141 
1 
1 
x=1+—=141 Z 5 
з 1+1 
1+ 1 
1 
1 8 
Xges T ld = etc. 
ё 1+1 
1+1 
1+1 
1 
It is interesting to see the Fibonacci numbers appearing as numerators and 
denominators of each x,, and we might make the conjecture that x, = Daa for 


n 
all positive integers n. This is certainly true when n = 1 and, if we assume that 


Uk + 
Xk = Ба then 
Uk 


1 1 Up; Uk _ Uk+2 
Хы = 1 + 1+ ; 
Хк Ик+1 Ик+1 Ик+1 
Uk 


n Я A u 
So mathematical induction confirms that it is true that x, = == for all n > 1. 


п 
The question remains as to whether or not the sequence {х„} produced by our 
iterative process converges. This question is the same as asking whether the 

ies 
sequence 77 
problem turns out to have a nice geometric interpretation. Imagine squares with 
edge-lengths 1, 1, 2, 3, 5, 8, 13,... . These can be joined, in order, to form a 
sequence of rectangles where, at each stage, the new rectangle is formed by 
adjoining the next square to the longer edge of the existing rectangle. This 
process is illustrated by the figure. Now the shape of each rectangle, as 
determined by the ratio of the longer to the shorter edge, is the ratio of 
successive Fibonacci numbers. So the question becomes: Does this sequence of 
rectangles have a limiting shape? 


of ratios of successive Fibonacci numbers converges, and this 


n 
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| E | 


8x8 


13x13 


5x5 


21 


The answer, as we shall see from the work of this unit, is yes. If, for the 
moment, we assume that the sequence has a limit k, say, then it is not too 
difficult to see what k must be. For sufficiently large rectangles the ratio of the 
lengths of the edges is very close to k, so we could take these lengths to be 
approximately ky and y. To obtain the next rectangle we add on a square of 
edge-length ky: 


= 
ky ky 
y ky 
= > 
ytky 


The ratio of the lengths of the edges of this rectangle is also approximately k, 
and so 


у + Ку 
ky 
Hence, if our sequence tends to a limit, it ought to be a root of the original 


=k or К?—К—1==0. 


{ à ios $ 1 5 
quadratic equation x? — x — 1 = 0. The positive root is k = + М5 a number 
1 2 


which is known as the golden ratio. It has the property that, given any rectangle 
with edge-lengths in this ratio, then adding a square to the longer edge 
produces another rectangle with edge-lengths in this same ratio. 


The golden ratio first arose from the problem of dividing a line segment AB at a 
: AB AC ,. | , 1 Sra 
point C such that AC “СВ this common ratio turning out to be Li Like 
another well known irrational number, л, the golden ratio continually crops up 
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in mathematics, and is particularly prominent in geometrical problems 
concerning regular plane and solid figures. As one example, it is the ratio of the 
radius of a circle to the edge-lengths of an inscribed regular decagon. There is 
little doubt that the golden ratio was used consciously by many Greek 
architects and sculptors, particularly in the structure of the Parthenon. The 
Greek theatre at Epidarros has 21 tiers of seats above the gangway and 34 
below, 21 and 34 being consecutive Fibonacci numbers. 

Our main purpose in giving this example is to illustrate what is meant by a 
continued fraction. Each of the expressions 


( 


1, 1+1, 1+1 , 1+1 ; ete 


is called a finite continued fraction. Their limit is an infinite continued fraction, 
and is written as 
1+1 
1+1 


This is the simplest and most well-known example of a continued fraction. In 
the more general infinite continued fraction 


а) +1 
a, +1 
a, +1 
dts 
the а; are positive integers, with the exception of ад which can be any integer. 


The theory of continued fractions is extensive and their area of applicability 


very wide, so in this unit we can hope only to give you something of the flavour 
of them. 


Section 7.1 is concerned with finite continued fractions and here we do much of 
the groundwork for subsequent development of properties of infinite continued 
fractions. Turning to infinite continued fractions in Section 7.2, the main 
problem confronting us is the question of convergence. We shall discover that 
every infinite continued fraction converges to an irrational number, and in 
Section 7.3 we settle the converse problem by showing that every irrational 
number can be expressed, in a unique way, as an infinite continued fraction. In 
Section 7.4, we break the irrational numbers down further by classifying, in 
terms of their continued fractions, those numbers which arise as irrational roots 
of quadratic equations. 


We are all in the habit of writing down, and even manipulating, expressions 


1+ ./5 , a Р 
such as л, log, 3, Lus etc., which represent irrational numbers. But when it 


comes to practical calculations we must deal with rational numbers. (Even 
today's computers are limited in how many decimal places of л they can 
handle.) So the question of finding good rational approximations to irrational 
numbers such as log, 3 is important. In Section 7.5, the final section of the unit, 
we see how the continued fraction of an irrational number provides us with 
‘best possible’ rational approximations. 


The unit is based on Sections 13.3 and 13.4 of B. 
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Finite Continued Fractions 


(i) 


(ii) 


(iii) 


READ B Section 13.3, page 299 as far as page 304, line 15. 


Page 300, line 11 

Because the numerators are not ones, Bombelli's example is not a 
continued fraction in the sense of Definition 13-1. It turns out that our 
insistence that the numerators must be ones does not impose any 
restriction, since every number has a continued fraction in this form. In 
Section 7.4 we shall find the ‘correct’ continued fraction for /13. 

It may be far from obvious to you why Bombelli’s example does give NAE 3. 


The following non-rigorous argument, which ignores all the real problems 


of convergence, might at least convince you that the result is correct. If we 
assume that 


х=3 +4 
6+4 
6+4 
6 +... 


then adding 3 to both sides gives 


x+3=6+4 
6+4 


But this is the denominator of the first 4 in the expression for x, so 
substituting we get 
4 
х=3+——— 
х+З 
which simplifies to x? = 13. Finally, since everything in sight is positive, we 
can discard the negative root and conclude x = AA 3; 


Page 301, line 3 

Any rational number, positive or negative, can be written with the 
denominator positive. Writing the number in this form before applying the 
Euclidean algorithm guarantees that only ao amongst the a; can possibly 
be negative. 


Page 302, line 7 
The complete Euclidean algorithm appropriate to finding the continued 
fraction for 19/51 should be 


19 — 0-51 + 19 
51 = 2-19 +13 
19— 1-134 6 
1322-641 
6-614 0. 


Notice that, like B, we have written the products on the right-hand side in 
our customary way for the Euclidean algorithm, but this reverses the order 
in which the products appear in the proof of Theorem 13-5. The 
underlined numbers are the a; of Theorem 13-5, and these are the partial 
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denominators of the continued fraction. Although B restricts the term 


‘partial denominator’ to a;,a5,..., we occasionally use it loosely to include 
do. SO p 
19 
s 9*1 
2+1 эре 
1+1 
2pm. 
6 А 
SAQ 1 B page 312, Problems 13.3, number 1. 
SAQ 2 B page 312, Problems 13.3, number 2. 
SAQ 3 B page 312, Problems 13.3, number 4. 
SOLUTION Using the notation of Note (iii) 
TO SAQ 1 
(а) -19 = —1-51 + 32 (b) 187 = 3:57 + 16 
51 21:32 + 19 5723-1649 
32— 1:19 + 13 16=1:9+7 
19= 1:13 +6 9=1:7+2 
13=2:6+1 7=3:2 +1 
6=6:1+0 2=2:1+0 
19 187 
So = [715 1, 1, 1, 2; 6]. So = [3:3, 151,3, 2]. 
(d) 118 = 0-303 + 118 
303 = 2:118 + 67 
(c) 71 2 1:55 + 16 118 = 1:67 + 51 
55 = 3۰16 +7 67 — 1:51 + 16 
16 = 2۰7 +2 51 = 3۰16 +3 
723:241 16 = 5:3 +1 
2=2۰1 +0 3=3:1+0 
71 118 
лыш [Г] А — = [0; 1 ; 
So 51 [173.2 3, 2] So 303 (0; 2, 1, 1, 3, 5, 3]. M 
SOLUTION (а) [—2; 2, 4, 6, 8] = –2 +I اک‎ 
TO SAQ 2 2+1 2+1 
4+1 4-48. 
6+1 49 
8 
204 710 
= —2 +1 =—2+— = -—, 
2345. 457 457 
204 
741 
4; 12. =—. 
(b) [4; 2, 1, 3, 1, 2, 4] 170 
321 
(c) [0; 1, 2, 3, 4, 3, 2, 1] = 460 


COMMENT You will no doubt be relieved to hear that we shall shortly 
discover a far simpler method of calculating these continued fractions. Bil 
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(a) (0; 3, 1, 2, 3] = [0; 3,1, 2,2,1), since 3=2+1. 
1 
(Б) [-15 2, 1, 6, 1] = [-152, 1, 7], зшсеб+1=7. 
1 
© 2;3,1,2,1,1,1]= [2;3,1,2,1,2] ѕіпсе1+1-2. NI 
1 


In SAQ 2 we computed continued fractions in what looks like the only logical 
way, starting with the final partial denominator and systematically working our 
way upward, calculating the various fractions along the way. In the next reading 
passage we learn an alternative method which starts at the other end, 
successively calculating the continued fractions in the sequence 


[ao], [403 a1], [do; ay, аз ],...,[ао; ay, 5, ..., ds]. 


As the continued fraction we want is the last term of this sequence, this may 
look like a lot of unnecessary work. It turns out, however, that each member of 
the sequence (from the third onwards) is easily calculated from the previous two 
without even involving fractions or common denominators. 


READ B from page 304, line — 10 to page 307, line 4. 


(i) Page 304, last line 
Although the formula given here is quite clearly valid, it must be noted 


that, except when a4, = 1, a, + 


is not an integer and consequently 
Oei 


ао; 04, d2,..., Qy 1, йк + is not a simple continued fraction. This 


A+ 1 
result is therefore set in the context of the more general continued fraction, 
where we adopt the same notation. The result is used later in the proof of 
the important Theorem 13-6. 


(ii) Page 306, line 12 
There is a slight flaw in the inductive proof of Theorem 13-6. In Note (i) 


we observed that |o: ai, d5,..., Am + is not in general a simple 


Om +1 


continued fraction, so the replacement of a,, by a,, + is not justified by 


т+1 
the inductive hypothesis. 
Fortunately, this flaw is easily remedied since the result of Theorem 13-6 
happens to be true for any continued fraction, simple or not. All we need 
do is omit the word ‘simple’ from the statement of the theorem and the 
words ‘the integers’ from line 8. The formulae for Co, С, and С, (on page 
305) which launch the induction, remain valid, and now the critical stage of 


the inductive step, replacing a,, by a, + , creates no problems. So we 


т+1 
have a valid proof of Theorem 13-6 for the more general continued 
fractions. As an immediate corollary, it holds for simple continued 
fractions, which are our immediate interest. A 


Notice that, for simple continued fractions, the denominator ак of the 
convergent C, is positive for К > 0. Certainly qo = 1 and а, = a, > 1 are 
both positive, and the inductive step is provided by the formula 

qk = 446-1 + Чк—2 for k = 2. ` 
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The formulae in the middle of page 305 are crucial. In what follows we shall use 
them many times, primarily for calculating convergents of a continued, fraction 
but also for theoretical considerations. The SAQs in the unit will provide you 
with plenty of opportunities to practise using them. 


In the following example we suggest a tabular method for calculating the 
successive numerators and denominators. This method may help you avoid 
arithmetic slips and might assist you in remembering these important formulae. 
After SAQ 4, which follows, we shall not use tables again, but you might like to 
persist with this method. 


Determine the rational number represented by 
[4; 2, 1, 3, 1, 2, 4]. 


In the tabular method we set out four rows labelled k, рр, a, and ак. The a, are 
the given partial denominators and so the entries for this row, and the k row, 
can be displayed immediately. The p, and q, rows are started by entering the 
first two values for each, using the formulae po = ao, py = aoa, + 1, qo = 1 and 
qı = ау. For our particular example these are po = 4, py = 9, qo = 1 and q, = 2 
and the initial table is 


> 


wl ea Ea tela 


а | 1 2 


To complete the table we move along the a, row filling in the entries above and 
below each a,. For the entry above a, (i.e. p,) we multiply a, by the entry 
diagonally above to the left, and add to this product the entry to the left of the 
latter. So, for example, p = (1:9) + 4 = 13, and then p, = (3:13) + 9 = 48. 


E. 


The same procedure applies for the q,. We multiply a, by the entry diagonally 
below to the left and add on the entry to the left of that. So q = (1:2) + 1 = 3, 
and so on. 
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The complete table is 


Re | LOR Ж une а AL 56 


Pe | 4 | 9 | 13 | 48 | 61 | 170} 741 


4|4|2|1|3]|1|2|4 


a |1]|2.3 | 11] 14] 39 170 


4 9 13 48 61 170 741 
and the successive convergents Pk are -, = d the last 


ge T2 ЗГ 39 Te 
being the value of the continued fraction. 


In our solution to SAQ 4 we have used the tabular method for parts (a) and 
(b). For part (c) we have used formula (*) of Theorem 13-6, calculating the 
successive numerators and denominators (and therefore the convergents) 
simultaneously. 


В page 312, Problems 13.3, number 5. 

(a) | E 
k|[0l1|2]|3/|41|5 
Px | 1 | 3 | 10 | 33 | 76 | 109 
ae | LA FSA al 
a | 1 | 2| 7 |23 | 53 | 76 
1.3 10 33 76 109 

The convergents are PP 7? 3353 and 76" 
(b) 


The convergents are 


(c) Using (*), the sequence of convergents of [0; 2, 4, 1, 8, 2] is 
01 (4:1) +0 4 (1:4) +1 5 (8:5) +4 44 
Г2 (4:2) +1 9 (1:9) +2 11 (8:11) +9 97 
(2:44) +5 _ 93 
(2-97) +11 205° 
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The next reading passage shows how continued fractions can be used to give a 
good method of solving linear Diophantine equations. d 


READ В from page 307, line 7 to page 310, line 15. 


B page 313, Problems 13.3, number 11. 


(a) The ane fraction for 35 is [2; 1, 2, 6] and the convergents are 2, 3, $ 
and 35. So by Theorem 13- 7, 51:3— 19.8 = (—1)27! = 1 giving one 
с а хо = — 8, уо = 3. The general solution is 


{х= —8 + 515 у 2 3 — 191: teZ}. 


COMMENT We could have started by finding the continued fraction for 
43, which would have been slightly longer since it has an extra partial 
denominator. Whichever method you use, make sure that you do not end 
up with the x and y values the wrong way around! 


(b) The continued fraction for 3$4 is [1; 1, 1, 1, 1, 10, 1, 3] and the convergents 
are 1, 2, 3, 5, 3, $3, 28 and 3$5. So 364-58 — 227-93 = 1, giving one solution 
хо = 58, yo = —93. The general solution is 


{x = 58 + 227t, y = —93 — 364t: teZ}. 


(c) We first consider the equation 18x + 5y = 1. The continued fraction for 28 is 
[3; 1, 1, 2] and the convergents are 3, +, 4 and 48. Therefore 
18:2 — 5:7 = 1, from which we obtain 18 - (2-24) — 5 : (7:24) = 24. This 
gives a particular solution xo = 48, yo = — 168 of our original equation, and 
consequently the general solution is 


{x = 48 + 5t, y = —168 — 18t: teZ}. 


(d) The continued fraction for 438 is [2; 1, 3, 2, 1, 1, 2] and the convergents are 
2, 1, M, 25 1$, & and 125. Therefore 158-22 — 57-61 = — 1. That is, 
158 :(— 22) – 57: (— 61) = 1 giving xy = — 22, yy = —61 as one solution. 
The general solution is 


{x = —22 + 57t, у= —61 + 158t:teZ}. Ш 
In the final reading passage of this section, we prove a result concerning the 
behaviour of the convergents of a continued fraction. This will play an 


important part in settling the question of convergence for infinite continued 
fractions, which we shall be turning to in the next section. 


READ В from page 310, line — 12 to the end of Section 13.3. 


Page 310, line —9 
The proof of the Lemma is not by induction at all. Having catered for the 
special case qo < 41, for k > 2 we have 


dk = Qkdk-1 + dk-2 > 6-1 (since q,.; > 0) 
and 
а4к-1 2 Чк-1 (since a, = 1). 


Therefore q, > q, ,,forallk22. A 


Calculate the convergents of the following continued fractions and in each case 
verify Theorem 13-8: 


(a) [13 2, 2, 2, 2, 2, 2,]; (b) [2; 2,4, 2, 4, 2, 4]. 
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(a) The convergents are Co = 1, С, 23,0, =1, С» 13, Са = $4, Cs = 2 and 
С, = 785. Theorem 13-8 says that Co < С, < © < C; «C5 «C4 «C, 
which we can confirm by evaluating the convergents to 4 decimal places: 

1 < 1.4 < 14138 < 1.4142 < 14143 < 14167 < 1.5. 


(b) The convergents are Cy = 2, С, = $, С. = ®, С, = $3, С, = 218 C, 133 
and Cg = 4438. Theorem 13-8 says that 
C, «C, «C, «C,« C, « C, « C,. 
In confirming this, to distinguish between C; and C, we have to approximate 
these rational numbers to 6 decimal places: 


2 < 2444444 < 2.449438 < 2.449489 < 2.449495 < 245 < 2.5. 


COMMENT A rather more elegant way of verifying that <$ is to calculate 


ad — be and show that it is positive. For the verification Cg < Cs we find 
485.881 — 2158-198 = 1, as we expect from Theorem 13-7. [ | 


Infinite Continued Fractions 


We now consider infinite continued fractions, for which we have a natural 
extension of the notation of the previous section, using (4o; 41, d2, a3,...] for 
the continued fraction with partial denominators (including ao) do, ау, аз, 
a3,.... We shall not use the adjective simple in connection with infinite 
continued fractions for it is to be understood that in any infinite continued 
fraction each of the a; is a positive integer, with the exception of ay which may 
be negative. Note, however, that the coming work necessarily involves finite 
continued fractions which are not simple. 
The first task confronting us concerns the question of convergence. Whereas the 
calculation of a finite continued fraction requires only the normal rules of 
arithmetic, the infinite continued fraction is an infinite process and as such 
requires a completely different approach. There is a strong resemblance to 

со 


infinite series. For the latter, we define the infinite sum У a, to be the limit of 
" r=1 
the sequence | Y a of partial sums, provided this limit exists, For infinite 
r=1 
continued fractions, the convergents play the róle of the partial sums and so, by 


analogy, we define the value of the infinite continued fraction to be the limit (if 
it exists) of the sequence {Cy} of its convergents. 


READ B Section 13.4, page 313 as far as the foot of page 315. 


(i) Page 314, line 12 
For the moment B shows only the existence of the limit, leaving the proof 
that it is irrational until later. 


(ii) Page 314, line —8 
By the Lemma, B page 310, the q,(k > 1) are a strictly increasing sequence 
of integers, and thus increase without bound. 


(ii) Page 315, line —7 
B also uses the notation [ao; ay, 45, ..., а„] for the infinite continued 
fraction [ao; а;,...,а, do, 05,. .., dn, do, ...]. (Notice that for this it is 
necessary that the integer ао be positive). Such periodic continued fractions 
are called purely periodic. A 
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The values of periodic fractions are readily calculated, as the following examples 
illustrate. 


We shall find the value of the purely periodic continued fraction 153,21], 


Та practice, we would proceed as follows. кеш x = [1;3,2] we observe that 
x = [1;3,2] = [1;3,2, 1, 3,2], so that 


х=1+1 
3+1 
2+1 


x 


This means that we can express x as the finite (though not simple) continued 


49 9x +4 
ti = uR ZN h hi 1,=,= 
fraction х = [1;3,2, х]. Now the convergents to this аге "7 and PEL and 
as the final convergent is equal to x itself, we have 
9x +4 
= р , ог 7x?— 6x 4-20. 
7x +3 


The positive root (which is clearly the one we want) of this quadratic is 


ыйы 


Тһе zi iatis for x in the first step of the above argument takes some 
liberties with properties of limits. To give the same argument correctly we 
would, instead, observe that 


C,=1+1 (for n > 3) 
341 
QE 
E 3 
which gives C, = Best tet Now from the definition of x = [1;3,2], 
1G. 4 3 


x = lim C, = lim C,_,, 


n> оо noc 


— : ; з 9x + 4 
so taking limits of both sides of the above equation gives x = AIT as 
before. : 


The underlying ideas of this example will work for any purely periodic 
continued fraction. For if x — [“0; 41, d2,..., a, ] (ag > 1), then 
X = [do; а, 45,..., а„ X]. With the usual notation, the final three convergents in 

п-1 Pn “Pn + Pn- 
р L р and хр Pn-1 
4n- ji "dn Xn + 4-1 
to x, we arrive at the quadratic q,x? + (qn-1 — p,)x — Pn—1 = 0 whose positive 
root is the value of the continued fraction. 

1 34437 2 
Consider y = [2; 5, 1,3,2] =2 + 541° where x = + as discovered in 
x 


this finite continued fraction are 


. As the latter is equal 


11 11 2 
Example 2. The convergents to y = [2;5,x] are 2, — eae and so 


5' SFE 
ut», 
_ x42 7 47 + 11,/37 


5x +1 3 + J17 22454531. 
5 T +1 
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Rationalizing, by multiplying top and bottom by 22 — 5, /37, this simplifies to 


y-18 - /37 


63 


This example illustrates how to calculate any periodic continued fraction 

Y = [405 41, а›,...,а,,, bi bo... b, ]. We first calculate x = [b, 3b2,b3,...,b,] by 
the method of Example 2 and then obtain у = [40; ау, a5,..., à,, x] by the 
method of Example 3. 


If you glance ahead to page 317 of B, you will find another example of this 
type, Example 13-4. 


SAQ 7 Evaluate the continued fractions: 


(а) [1:2], (b) [1:12] © [2;1,1,1,4], (d) [—2;1,2,3,4]. 


OLUTION (0 Lei (03) — eT apd ga 2991 
aos YE 2x +1 
X 
2x* — 2x — 1 — 0, which has positive root x — ом 


Therefore 


= 1 1 
(b) Let y = [1;1,2] = Ш;х]=1+ v where x — 8 from (a). Therefore 
— 343 1-3 Ji 
CORBIS 14/5 1-8 
(c) Let x = [1;1,1,4] = [1;1,1,4, x]. The convergents are 1, = 5 5 and 


1 
14x + 3 14x +3 
м = ЕЕ 2 کا‎ x 
Ge 15 herefore x 9x 42” giving 3x^ — 4x — 1 — 0, whose positive 


root is x — mo 


Now, 


——— 1 3 
(231,1,1,4] = [2;х] 2 24 =2+ E A T. 
1= [23x] = nA JI 


13x +3 


= [3:4] = [3: =. 1 =T р 
(d) Let x = [3;4] [3;4,х] 341 Sox ах 1 Bing 


Now, 


LI 4x-1 -14-8,/3  -10-4,/3 
[72:,2,3,4] = [-2:,2, x] = Шш 
3x +1 11 +63 13 


In the final reading passage of this section we prove two important results 
concerning infinite continued fractions. First we confirm that the value of any 
infinite continued fraction is an irrational number; and second, in anticipation 
of the converse result, we prove that the continued fraction of an irrational 
number is unique. 


READ B from the top of page 316 as far as page 317, line — 13. 
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Page 316, bottom 
The observation made here can be extended in an obvious way: 


[40;41,45,...] = ао +1 
[413 d2, 45, ...] 


=a +1 
a,+1 = 
[42; 45, d4,...] 
and these results provide a rigorous basis for the intuitive arguments used in 


our Examples 2 and 3. B is quite right to give justification, for when dealing 
with infinite sums, nothing is obvious! A 


B page 327, Problems 13.4, number 2. 


Consider [0; ao, a1, a5,...]. (Notice that x > 1 ensures that ao = 1 and so this is 
an infinite continued fraction.) By the result at the foot of page 316, 


[0; ao, а„,а›,...]=0+1 Lu 
[40;41, 42, ...] 


5 ; = — 1 
Using this result with x = [1;1] we see that [0;1] = ту However, we have 
= 1 5 
already seen on page 315 that [1;1] = Lint so 


[0;1] 2 y5-1 п 


ETT EE 


The Continued Fraction Algorithm 


We have now seen that every infinite continued fraction represents an irrational 
number. We are left with the question of whether every irrational number has 
an infinite continued fraction. In this section we show that the answer is yes, 
and Theorem 13-10 confirms that the infinite continued fraction of an irrational 
number is unique. 


If, for the moment, we assume that the irrational number хо has an infinite 
continued fraction, say xo = [ao; d, а, аз,...,}, then the proof of 
Theorem 13-10 tells us exactly what the a; must be. 


We know already that aj = [хо], the integer part of xo, and furthermore 
1 1 


Xo = do + [xo] + E 
4 s [41;42, 45,...] ? [41 42, 5, ...] 
1 
But this gives ————~ = [a,; a, d3,...]. 
8 xo — [xo] [415 42, ds, ...] 
Р 1 ; К AM TEN 
Now if we put x, = —— which, like хо, is irrational, then 
Xo — [xo] 


X, = [a,;425,45,...]. As before, a, = [x,] and 


1 
х= а і [xı] + 
! á [a2;a3,...] ^ [а2;аз,...] 
; 1 — ; 
And so we can continue, x; = Ec (an irrational), which leads to 
1— 1х1 


аз = [x5], and so on. 


56 


EXAMPLE 4 


M335 7.3 


In fact we have a procedure (known as the Continued Fraction Algorithm) for 
computing successive partial denominators for the continued fraction of Xo, 
which is summed up in the following flow chart. 


Irrational x, 


Express this number as integer part + remainder 


Calculate the reciprocal 

partial denominator of the remainder 
Of course, it still remains to show that the infinite continued fraction 
determined by this process converges to the original irrational number хо. The 
proof of this is contained in the next reading passage from B, but before that let 
us assume the result and illustrate the algorithm in action. 


In Section 7.1, we saw Bombelli’s continued fraction for ЫЙ. We shall now find 
the ‘correct’ one. 


Suppose 4/13 has continued fraction [ao;a;, а, аз,...]. We follow the above 


algorithm. 
хо = 4/13 = 3.6055..., 
so do = [Xo] = 3 and xo — [xo] = 4/13 — 3. 
1 1 1 f/13 +3 /13 +3 
Xi . > 
Xo — [xo] 13-3 13-3 \/13 +3 4 
6.6055... _ | 13 =1 
4 4 p 
J13 —1 
soa, = 1 and x, — [xi] = Т 
1 4 J13+1 4.6055... 1+ [=] 
x T > 
* Ha nl fo <4 3 3 3 
13-2 
so а, = land x, — [х] = 3 à 
1 3 J13 +2 | SS 
x E 
° Xa De] {13-2 3 3 
J13—1 
so аз = land x4 — [x3] = 3 Я 
1 3 /13 +1 MERE 
X4 1 i 
хз — [x3] 13-1 4 4 
J13—3 
so аа = 1 and x4 — [x4] = 4 ; 
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1 4 
v13 +3 = 6.6055... = 6 + (,/13 — 3), 
ха — [X4] 4/18 —3 ыз ) 


so as = 6 and х; — [xs] = 4/13 — 3. 


1 " 1 
—= = x, implying that a = ау. Moreover; x, = — — _ 
/B=3 1 plying 6 1 ym xe — [xe] 

1 


= aoe = x3 gives a; = az. Clearly, we are in a cycle, and the sequence of 
Xi — [Xi 


Xs 


But now x; — 


partial denominators 1, 1, 1, 1, 6 will recur indefinitely. So, assuming that ./13 
can be expressed as an infinite continued fraction, it must be [3; 1, 1,1, 1,6 i^ 


READB from page 318, line 5 to page 323, line 9. 


(i) Page 319, line 11 
Recall the convention for convergents. C, is called the nth convergent. 
However, as C, is the zeroth convergent, there are in all, up to and 
including C,, n + 1 convergents. 
On line 14, C4, , is the last convergent of the finite continued fraction on 
the previous line, and as such is the value of the fraction itself. 


(ii) Page 319, line — 10 


Once again, the argument here requires our stronger form of the fraction itself. 


Theorem 13-6, where the contined fraction is not required to be simple. А 


B page 327, Problems 13.4, number 4 parts (a), (d) and (e). 
B page 327, Problems 13.4, number 5. 


(а) \/5 = 2236... 
хо = /5 =2 + (5 — 2), do —2 
ma Tea ل‎ +2 =4 + 5 2), а=4 


E 
9 


Therefore „05 = [234]. 


(d) 4/37 = 6.083... 


= xı, and so we are in a cycle with 4 = a; = a, =a; =... . 


5 / = 
Xo = = E 3! ag = 2 
4 4 
4 J31 +3 J[31 —4 
Xi Tp š а; =1 
37-3 7 7 
7 37 +4 J31—5 
X? 3+ b a, =3 
37-4 3 3 
Xz = = 2 Y 1+5 Xo, and so we are іп a cycle starting with 


аз = ао. Therefore tas 37 = [2;1,3]. 
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\/30 = 5471... 
„Шын. 1+ (4), "em 
13 
+2 т = 
1=3 
“Tina 2 
" 2 „30 + 4 r+ م‎ m 
2 30-4 7 2 
— | B0 +3 2+ (09-3) m 
P 3 3 
xd _ /3043 ie p, -— 
353 3 7 7 | 
—: 30 +4 p) xix d 
/30 — 4 2 2 | 


Xe = —у=——— = X» SO we are in a cycle starting with ag = ay. 
/30 — 4 | 


Therefore иын = [13,1214]. M 


Using the hint, we see that /n? +1 =n + 


1 
— 7———— where, clearly, 
n + Jn? +1 


0< a < 1. Therefore 
n + /n*+1 


1 
ж= +1 =n [LÀ 7 а =n 


xı =n + Jn? +1 = 2n + (n? +1 - п) = ar TIT a, = 2n 


x =n + Jn? +1 = х, and so we are in a cycle with 2n = a, = d= 


Hence /n +1 = [n;2n]. 


xo = Jn + =n + ((/n? +2 — n), where 
0< Jn? +2 —n ae | do =n 
n 


2+2 +n 


a =n 


1 n? +2 +n ا‎ 
x= n+ 2 > 


Jn! +2=n 2 


ee ee eee Fu NE a, = 2n 
Jn? +2 —n 

i Xı, so we are in a cycle starting with аз = a 
хз =—= = x; t = a. 


Hence \/n? + 2 = [n;n2n]. 
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Мп + 2n = n + (mh + 2n — n), where 


veu 2 
0< M ger — 
Jn? + 2n +n 


Xo =n + (\/n? + 2n — п), do = n 


| FER A 
x —— 
t Jr + 2n п 2n 
-14 Vta 
2n 
1 
=1+-=———— a,=1 
Jn + 2n +n 
x, = Jn? + 2n + n = 2n + (Jn + 2n — n), а, = 2n 
1 


Хз = Ta = Xa, SO we are in a cycle starting with аз = a,. Hence 
Jn? + 2п —n 


n? + 2n = [n;l2n]. 
(b) /2= JP? +1 = [152] 


3= /1 +2 = [1512] 
15 = 4/32 + 2۰3 = [3;1,6] 


37 = /6 +1 = [6;12] M 


The Corollary оп В page 320 tells us that, with the usual notation for the 
convergents of the irrational number x, 


1 
annsi 


This can be used to predict how accurately the convergents approximate x, or 
alternatively to select the appropriate convergent for a desired degree of 


3.22. 333 355 
accuracy. We have seen that the first five convergents of x are T 7° 106 13 
103993 2 
and 331027 From the Corollary, |x — 7 < 7-106 ^ 0.0013... < 0.005, so 


certainly 3 gives л accurately to two decimal places. On the other hand, if we 
want an approximation to z which is accurate to six decimal places, we have 


« 0.0000005 to see that эз is such an 


1 
only to observe that 11333102 133 


approximation. 


SAQ 11 В page 327, Problems 13.4, number 6. 
SAQ 12 B page 328, Problems 13.4, number 7. 


SOLUTION (/15 = [3;1,6] has convergents 


3 
TO SAQ 11 1 


4 27 31 213 244 1677 
1027 8 NSS 63? 43377" 


> 


1 244 
= = 0.0 6... < 0. 5 —i four decimal 
6-433 0.000036... < 0.00005, the convergent a 5 correct to four decima 
places. B 
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(a) The convergents are 


., and as 


2 3 8 11 19 87 106 193 1264 

113 4° 7°32 39° 71" 465°” 
193 . 

"Id 0.000030... < 0.00005, the convergent up correct to four 


decimal places. 


a 87 
(b) If e *y < a» then 


e 


87 1 
32| © 32-39 
. 106. | T 

by the Corollary, since 39 is the convergent following a So 
87b — 32a 1 
— < 

32b 32:39 
and as 87b — 32a > 1, we must have b > 39, IB 


0< 


The Continued Fractions of Quadratic Irrationals 


The quadratic irrationals are the numbers of the form r + s / d, where r and s 
are rational, s # 0 and d > 0 is a non-square integer. They are the numbers 
which arise as irrational roots of quadratic equations with integer coefficients. 
On page 322, B remarks that a number x is a quadratic irrational if and only if 
the continued fraction of x becomes periodic ultimately. In this section, we shall 
examine this claim, proving the ‘if’ part. We begin, however, by considering the 
simpler case of the continued fractions of Jad, where d is a non-square integer. 


At the bottom of page 334 of B, there is a table giving the continued fraction of 


Jad for every non-square integer in the range 1 < d < 40. Concealing the 
remainder of page 334, examine the given table and then have a go at answering 
the following. 


What common features are there in the continued fractions given in the table on 


page 334? Conjecture the general form for the continued fraction of Ja, where 
d is a non-square positive integer. 


From the evidence presented you might have made some of the following 
observations. 


(i) The continued fraction of Jd is periodic, with period starting at a, 

(ii) The final element in the period is 2ag. 

(ii) Leaving aside the final term, 2ao, the remainder of the period is symmetric 
(palindromic), reading the same forwards as backwards. In some cases, 
such as Ja, „5 and 4/10, this symmetric part is non-existent. 


Collecting all this together, we might conjecture that the continued fraction of 
Jd takes the form 


d = [а0; dı, а, аз,..., аз, аз, ау, 2d0]. 


COMMENT We can express this conjecture in an equivalent and more 
succinct form by observing that if it is true, then [v4] = ао so that 


[V/d] + „= [240: %1, 45, 5... ., A3, 5, A}, 245 ] 


= [2a0;a1, 49, 43,...,43, 42,41 ] 
which is purely periodic. Bil 
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By the end of this section we shall have gone some way towards proving the 
conjecture in SAQ 13. As a start in that direction, we investigate the 
relationship between a purely periodic continued fraction and the fraction 
obtained by reversing the order of the elements in its cycle. 


Let « = [2;1,4,1] and f = [1;4,1,2] (obtained by reversing the period of a). 
How, if at all, are х and В related? 


= 2-8 
Writing « = [2;1, 4,1] = [2;1, 4, 1, а], the convergents of the latter are TT 
14 17 17a + 14 , БЕ 
56 and E æ, SO that « is the positive root of 
6x? — 120 — 14 = 0. (1) 
——— 6 17 17 M 
Similarly, f = [1;4, 1,2 ] has convergents 1, E $ : and ju a В, giving 
148 — 128 -6=0 (2) 


of which f is the positive root. 
One cannot help noticing the similarity between the coefficients of the 


quadratics (1) and (2). Indeed, if we substitute x = Zd in (1) we get (2). This 


В 
—1 
means that « and d are the two roots of (1). (They cannot be the same root 
since « > 0 and F < 0.) We can check that this is so by actually solving the 


: 3+ 30 з o ; 
equations. The roots of (1) are =a being the positive one, whilst the 


positive root of (2) is В => У 30 Now S 37 30 


B 344/30 3 


The pair of quadratic irrationals r + s/d {which are the two roots of some 
quadratic equation with integer coefficients) are called algebraic conjugates. So а 


„аз 


claimed. 


' ==, . : 
and f are related, in that E is the algebraic conjugate of x. 


Let x = [2;3,5,1] and $ = [1;5,3,2]. Show that о and y are algebraic 


conjugates. 
2 7 37 44 44a + 37 
The convergents of x are T, 16 19 19:416 ^ and so 
19a? — 280 — 37 = 0. (1) 


161944 E d 
CS 16 37-376 4-16. ^» 20080 


374? — 288 — 19 — 0. Q) 


The convergents of f are 


meh. cu А : -1 "e 
Substituting E for о in (1) gives (2), and so х and 7 are the two (irrational) 


roots of (1), and as such are algebraic conjugates. 


—1 . 144 8 А 
(For the record, x and d are the pair ae but notice that there was no 


need for us to calculate х and f.) Wi 
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In Example 5 and SAQ 14, the all-important relationship between the 
coefficients of the quadratic equations in x and В stemmed from a connection 
between the last pair of convergents of x and fl. In Example 5, we had 

14 17 6 17 37 44 19 44 
Eur for « and =, 14 for B. In SAQ 14, we had 16 19 for « and ТЄ 37 for f. 
This suggests the following general result, which we prove on our way to a 
classification of those irrational numbers which have a purely periodic 
continued fraction. 


Let the final two convergents of [а0;а;,а›,...,а,| be Pat and Pa. Then the 
q 


n-i n 


Jw. and Pa, 
@һ-1 Pn-1 


final two convergents of [a,;a,..,...,a,,a9] are 


(By mathematical induction on n.) 


А Р . a 

The case п = 1 is easily verified. The only two convergents of [a9;a, ] are Y 

dod; +1 > а ава, +1 

кесе, whilst the two convergents of [a,;a9] are 7 and =" The 
1 do 


pairs are related as claimed. 


and 


The inductive step is a little more complicated, but provides good opportunity 
to practise manipulation of convergents. 


Complete the proof of the Lemma. 


Assume the result is true for n — k. With the usual notation for convergents this 

means that the final two convergents of [ao; 41, d2, ..., ap] are = and 5 and 
к= 1 [4 

Че. апа Ры. 

dk-1 Px-1 


the final two convergents of [a,;a,..,,...,d,,a9] are 


> р 
Consider [40: @1,..., 454], whose final two convergents аге = and 
k 


a » 
Prva _ Geri Pu t Рк-1 We need to prove that the final two convergents of 


Яка +10 + 9-1 


йк+1Чк + Чк-\ and ак+1рк 


+ Py Р 
[дөге йу] are Рк-1 That is, we must 


4, Pk 
show that 
А __ 419k + 9-1 
[ак+1;аь...,ау] = ——— — — 
qk 
and 
E ак+1рк + Px-1 
[0+ 1; +++, 4, 49] = —— —— ——— 
Pk 
Now 


1 
[ay 14... 2,41] = Gi + 
лан н [ак;ак-1,...›@2,ау] 


1 ak+14k T Чк-1 
d dk -1 dx 


(Note that [а,; а, – 1,...,а,а; ] = 8/4: by the inductive hypothesis, being the 
penultimate convergent of [a,;a,—1,.-.,41, 4o ].) And 


йк+ + 


1 
Ay 1; Ay. ++, 01, а+1 + 
[@k 151 +++, 41, 00] En Tae ас 


1 ак+1рк + рк-1 ш 
рьрь-1 Рк 


Ars + 


63 


THEOREM 1 


CONVERSE 1 


THEOREM 2 
PROOF 


CONVERSE 2 


M335 7.4 


Returning now to our purely periodic continued fraction, we shall establish the 
general result which we illustrated in Example 5 and SAQ 14. 


If «= E TE ], then (see Example 2) о is the positive root of the 
quadratic equation 4,0? + (q,_; — p,)x — р: = 0. Reversing the order of the 
elements in the cycle of 2, we get 

B= [а,;а,-1,...,а1,а0] = [аһ;а„- 1,...,d1, 40, В]. By the Lemma, the final three 


4н 3 P" and Bp, + d, B. Hence f is the 
Qn-1 Pn-1 BPn-1 + 4-1 


convergents of the latter аге 


; e]. 
positive root of р, 12 + (dn—1 — Pn) B — q, = 0. Replacing x by F in the first 
—1 У ; 
of these quadratics gives the second, and so « and E are algebraic conjugates. 


+ 1 
Collecting all this together, noting that x > 1 (since ао > 1) and -7 > —1, we 


have proved: 


A purely periodic continued fraction is a quadratic irrational greater than 1 
whose algebraic conjugate lies between — 1 and 0. 


It turns out that the converse of Theorem 1 is also true, but this is a much 
deeper result, which we are not going to prove here. The first proof was given 
by Galois in 1828 although, from earlier works, the result seems to have been 
known to Lagrange. 


Any quadratic irrational greater than 1 whose algebraic conjugate lies between 
—1 and 0 has a purely periodic continued fraction. 


As a consequence of Theorem 1 and its converse, we can also classify all 
periodic continued fractions. 


Any periodic continued fraction is a quadratic irrational. 

Let y = [do341,...,4n,b0,b1,...,b, ] be periodic. Then x = [bo;b,,...,b,,] is 
purely periodic and so, by Theorem 1, x is the positive root of a quadratic 
rx? + sx + t = 0, where r, s and t are integers. 


Writing y = [a0;4;,...,a,, X), with the usual notation for convergents 


nX + Pu-i Q2. n-1Y — Pn- 
m р L giving x = 4 1Y — Pn-1 
@һХ + qn-1 Pn — day 
quadratic we obtain 


. Substituting this value of x into the 


P 


y) — 2 f 4 
4-15 а) ere bia 
Pu — 4У Pn — у 


| "E 
or 
F(Qn-1Y — Pa-1)” + 5(4„-1У — Pn—1) (Pn — qay)  t(p, — Guy)? = 0 


which is of the form Ау? + By + C = 0, where A, В and C are integers. As y 
must be irrational, this confirms that y is a quadratic irrational, as claimed. $ 


The converse of Theorem 2 is also true. Its proof follows without too much 
difficulty from Converse 1, but again we shall omit it. 


Any quadratic irrational has a periodic continued fraction. 


In conclusion, let us apply these results to Jd, where d is a non-square integer, 
to verify our conjecture of SAQ 13. 


[Jd J+ Jd is a quadratic irrational greater than 1 whose algebraic conjugate 
[V4] — \/d satisfies —1 < [V/d] — \/d < 0. Hence, according to Converse 1, 
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[V4] + Jd has a purely periodic continued fraction, say 
[V4] F Jd = [240;4,,45,...,4,], 


where dy = L/dl, since the integer part of L/d] + Jd is 2 [V4]. Reversing the 
period gives a continued fraction whose value is minus the reciprocal of the 


algebraic conjugate of [V4] + afd. That is, 


1 1 
[4n3 Qn— 15+ +++ 41,249] FE (1) 
i : Ë o Jd Jd — do 
But 
i 1 
Ja + do = 2a) + (va ao) = 2ay + - 1 
La - z 
and 
n —— шшш: 1 
Vd + ао = [2а0;а1,...,а,] = 2а, + — ——— 
[t Tbe, An, 200] 
Therefore 
1 ———— 
= = [a1;45,..., Am 200]. (2) 
Ja — ao 
By the uniqueness of continued fractions, comparison of the expressions (1) and 
1 ; ad , 
(2) for —— gives d, = а„ d = Ap- |, etc., establishing the symmetric 
[c 2 n g 
V TUM 


property of the period. 
Collating this information. 
[V4] + Jd = [2ао;а\,а»,.... аа) 
ог, equivalently, 
МЗ = [ao;41.45,..., 4, 41,200 J. 
Which of the following quadratic irrationals have a purely periodic continued 


fraction? 


2 7 1 37 
(a) pel, (itx q 5.24. 


6 


2464/7 et | А 
(а) TERT does not have a purely periodic continued fraction because its 


— 6/7 
5 


2 
algebraic conjugate, — 


1 37 
HET 


has a purely periodic continued fraction. 


, is less than — 1. 


. " А 1 — ./37 i 
> 1 and its algebraic conjugate satisfies —1 < EET NES « 0, so it 


(c) Although the algebraic conjugate —7 4- 2/11 lies between —1 and 0, since 


—7—2 Jil is not greater than 1, it does not have a purely periodic 
continued fraction. B 


Rational Approximations 


So far we have seen that any irrational number has a unique infinite continued 
fraction and that we obtain a sequence of rational approximations to the 
irrational number by taking the convergents. We even have an estimate for how 
good these approximations are. However, we have not yet fulfilled the promise 
made in the introduction that the continued fraction provides what are, in some 
sense, best possible rational approximations. 
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Given an irrational number x, we can always improve on any rational 
a 


approximation 5 In fact, there are infinitely many rational numbers lying 


a 


between b and x. Consequently there is no such thing as the best rational 


approximation to x. However, if we impose restrictions on the size of the 
denominator, we can ask about the best rational approximation to x. For 
example, of all rational numbers whose (positive) denominator does not exceed 
50, the best approximation to z is 22. (To improve on 22 we have to go to 7.) 
We shall show that if p,/q, is any convergent of x, then there is no rational 
number a/b with 0 < b < q, which approximates x better than р/а, and in this 
sense the convergents of x are best rational approximations to it. 


Before embarking on the reading passage we must give you a beautiful 
geometric interpretation of the continued fraction for an irrational number, 
which was first given by Klein in 1895. It illustrates precisely the sense in which 
convergents are best rational approximations. 


Let о be a positive irrational number and consider the part of the line y = ax in 
the first quadrant of the Cartesian plane. This line cannot pass through any of 
the lattice points (i.e. points whose coordinates are both integers) except (0,0), 
for if it did, x would be rational. Imagine that pins are set upright at each of the 
lattice points in the first quadrant, and a fine string is placed along the line with 
one end fixed to the line at an infinitely remote point, and the other end at the 
origin. Keeping the string tight we move the end at the origin towards the right. 
It will catch on a series of pins below the line. In fact, the pins on which it 
catches are those at lattice points with coordinates (do. Po). (d2. P2). (Q4. P4). ... 
corresponding to the even convergents po/qo, P2/42, P4/44,... Of a. 


Similarly, if we move the string from the origin to the left it catches on pegs at 
(d1 P1). (d3. P3). (5. P5). ... corresponding to the odd convergents p;/q1, p3/d5, 
P5/45,... Of a. 

813 


134 . ° . ° . ° 


, 1 5 Р 
The figure shows this for the case х = ке the golden ratio, whose 
convergents are f, 7, 3, 3, $, ,... . The dashed line is the portion of у= ax in 
the first quadrant. Moving to the right it catches on the pins with coordinates 
(1, 1), (2,3), (5, 8),... whilst moving left it catches on the pins at (1,2), (3,5), 
(8,13... 
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So much for the preview. The proof of the first Lemma in the reading passage is 
somewhat involved, but once that is out of the way the rest is plain sailing. 


| READ B from page 323, line 10 to page 325, line 10. | 


SAQ 17 (a) Prove that the convergents of x (irrational) are successively closer to x. That 
is. 


"m b E 
(n1 dn 
(b) Prove that the convergents of x satisfy 
1 1 
& x ا‎ ә . 
24,4+1 d 4+1 


SOLUTION (а) This is an immediate consequence of Theorem 13-12. Since the q, are 
TO SAQ 17 increasing, q, < q,., and so 


x= Bett اک‎ Pu 
Qni dn 
Equality cannot occur, for if it did, we would have 
х= Pia Be x (since x lies between successive convergents) 
qn+1 Qn 
; _ 1 [Pn _ Pasi ا‎ WOES 
which leads to x ERE ; contradicting the irrationality of x. 
Qn an+ 


(b) The right-hand part of the inequality is the Corollary on page 320 of B. For 
the left-hand part, note that 


1 Pn ра+1 К Pal , : Disi 
Qn 1 dn An+1 Qn Яп+1 
(віпсе[р,д, 41 — Р„+1йһ| = 1 and x — +", x — Pr*1 haye opposite signs). 
n ntl 
Now, using part (a), 
«2x — P» 
Gun 1 dn 


which gives the required inequality. Ш 


The final reading passage provides us with a test for whether a rational 
approximation to x will occur as one of its convergents. 


READ B from page 325, line 11 to the end of Section 13.4. 


SAQ 18 В page 328, Problems 13.4, number 10. 


SOLUTION 1 + /10 18 _ 13/10-41 (since 132-10 — 41? = 9) 
TO SAQ 18 3 13 39 13(13,/10 + 41) 


18 3 P 1 
13 13:80 2-134 


Now 4/10 3, and so : =» = 


18. 1 10 

Theorem 13-13 establishes that n is a convergent of LESER 

1-4 10 
3 


1347 18 
1,2,3 513- BH 


For verification, has continued fraction [1;2, 1, 1 ] which has convergents 
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Summary 


In this unit we have studied continued fractions and some of their properties: We 
started by investigating finite continued fractions, many of whose properties 
were important for the subsequent development of the infinite continued 
fraction. We saw how to use the Division Algorithm to obtain a continued 
fraction for any rational number, and a generalization of this technique, the 
Continued Fraction Algorithm, led us to the unique infinite continued fraction 
of an irrational number. 


One of the themes of the unit was to classify the numbers associated with 
various kinds of continued fraction. To this end we saw that rational numbers 
correspond to finite continued fractions, irrational numbers to infinite continued 
fractions and quadratic irrationals to periodic continued fractions. The latter 
category was broken down further by classifying which quadratic irrationals 
correspond to the purely periodic continued fractions. We also found the 
particular form of the continued fraction of Jd, where d is any non-square 
integer. 


We discovered a good method, based on continued fractions, for solving linear 
Diophantine equations. (In Unit 8 we shall use continued fractions to solve 
another family of Diophantine equations.) 


One of the main applications of continued fractions is in the area of finding 
rational approximations to irrationals. We discovered that the convergents of an 
irrational number are, in a certain sense, best rational approximations. 


The main definitions and results encountered in this unit are listed in the 
corresponding entry in the Handbook. 
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Introduction 


Our first encounter with a Diophantine equation was in Unit 1, where we were 
confronted with the linear Diophantine equation ax + by + с = 0. Prior to this 
course, you would have probably associated with this equation the geometric 
picture of a straight line in the plane, and the solution set of the equation 
consisting of those ordered pairs of real numbers (x, y) corresponding to points 
on the line. But in number theory we restrict ourselves to the integers, rather 
than real numbers, so the solutions we require are the ordered pairs of integers 
(x,y) which satisfy the equation. Geometrically, we seek the lattice points which 
the line ax + by + с = 0 passes through. There are a number of questions that 
we can ask, and have asked, of the linear Diophantine equation: 

(i) Are there any integer solutions? 

(ii) Is the set of integer solutions infinite? 

(iii) Is there a method of finding all integer solutions? 

We had little difficulty supplying answers to all three questions for the linear 
Diophantine equation. But the same questions can be asked of any polynomial 
equation P(x;,x5,...,x,) = 0 inn variables, and, more often than not, these 
questions present considerable difficulty, even when n = 2 and P isa fairly 
simple polynomial. For example, consider x? — y? — 2. You might be able to 
answer question (i) by spotting a solution, but how about questions (ii) and 
(iii)? (Fermat managed it, but we shall not go into his solution here.) 


Questions of this nature are called Diophantine problems, named after 
Diophantus of Alexandria, who initiated their study through several problems 
involving squares and cubes which arose out of geometrical considerations. 


Number theory is a vast subject, and there are many branches which we have 
not touched upon in this course. However, nearly all branches have as their 
motivation the objective of making a contribution to the solution of one or 
more Diophantine equations. It is this topic of Diophantine equations which, 
historically, lies at the heart of number theory, and the study of which has, over 
the ages, stimulated to a considerable extent the development of mathematical 
thought. 


In this unit we are going to look at a handful of the many notorious 
Diophantine equations, each one of which, as you will read, is surrounded by its 
own fascinating history. The material is selected from B, Chapters 11, 12 and 13. 


In Section 8.1 we continue where we left Unit 7, looking at an application of 
continued fractions. We shall find that Pell’s Equation, x2 — d y? = 1, where the 
positive integer d is not a square, has infinitely many solutions, all of which are 
readily obtained from the continued fraction of Jd. 


Section 8.2 investigates the Pythagorean equations. Perhaps the best known 
theorem in the whole of mathematics is that of Pythagoras, which says that if x, 
y and z are the lengths of the sides of a right-angled triangle, z being the length 
of the hypotenuse, then х2 + y? = 22, Solving this as a Diophantine equation 
amounts to finding all right-angled triangles whose side lengths are integers. 


Section 8.3 turns to the equations x" + y" = 2", for n> 2. These equations hold 
a prominent place in mathematics because of the history surrounding 
unsuccessful attempts to show that they have no solutions (other than the trivial 
ones in which x or y is equal to 0). The conjecture that they have no solutions 
has become known as Fermat's Last Theorem. We shall use our solution to the 
Pythagorean equation to prove Fermat's Last Theorem for the special case 

n — 4. This proof involves us in the one new technique of this unit: Fermat's 
method of infinite descent. 


Finally, Section 8.4 looks at the problem of representing integers as sums of 
squares. We shall give a complete solution of the two squares problem, that is, 
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we shall classify those integers N for which the Diophantine equation 

x” + y = N has solutions. The section rounds off with a brief look at the allied 
three squares and four squares problems and includes the result that every 
positive integer is a sum of four squares. 


Pell’s Equation 


We start our survey of Diophantine equations with one which has a curious 
history, Pell’s equation x? — dy? = 1, where the positive integer d is a non- 
square. 


2 
As Jd is irrational, the equation — = d, ie. x? — d y? = 0, has no integer 
y 


2 
У PE : — 91а 
solutions. But if — is а good rational approximation to Jad, then — will be 


relatively close to d, or what amounts to the same thing, x? — dy? = k will be a 
small integer. As k cannot equal 0, the next best thing is that k = +1. Having 
seen that the convergents of the continued fraction of Jd are ‘best’ rational 
approximations, it will come as no surprise to discover that every solution (х,у) 


; А ; x 
of either of the equations x? — dy? = +1 arises as a convergent — of Jad. 
F 


The bulk of the first reading passage is concerned with the origins of Pell’s 
equation. 


READ B Section 13.5, page 329 as far as the bottom of page 333. 


B page 341, Problems 13.5, number 3. 


Knowing that all solutions are to be found amongst the convergents of the 
continued fraction of m d, we need check, in each case, only those convergents 
with denominator less than 250. 

(а) x? – 2у2=1 


The continued fraction for КА is [1;2] and so the convergents are +, 3, 1, 
12 29 46> tes... . All subsequent ones have denominator exceeding 250. The 
corresponding values of pf — 2q; are, respectively, —1, 1, — 1, 1, — 1, 1, 
—1,..., and so the required solutions are (3, 2), (17, 12) and (99, 70). 
COMMENT Note how we can use Theorem 13-15 to assist in calculating 
p; — 242. It tells us that |p? — 242] < 1 +2 V/2 < 4, and so we need only look 
at the units digit of p; — 242. For example, for the convergent 229, 


239? — 201692 29? 2.92= 12.1 = 1 = 9 (mod 10) 
and so we must have 239? — 2-169? = —1. 


(b) x? - 3y? =1 


The continued fraction for WE is [1;1,2], and the convergents with 
denominator q, < 250 are 1, 7, $, 2, 19, 2% 21, 22, 265 and 382. The 
corresponding values of р; — 342 are respectively —2, 1, —2, 1, —2, 1, —2, 
1, —2 and 1, and the required solutions are (2,1), (7, 4), (26, 15), (97, 56) and 
(362, 209). 


(с) x? — 5у2 = 1 


From the continued fraction J5 = [2:4], the convergents of 5 are 2, 2, 38, 
3, ..., with corresponding values of p? — 542 being —1, 1, —1, 1,... . The 


required solutions аге (9,4) and (161,72). B 
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We 


now know that any solution of Pell’s equation x? — dy? = 1 must arise as 


the numerator and denominator of a convergent of Jd. What is more, in each 

of the examples, we have seen that solutions always occurred. Will this be true 

of every non-square positive integer d, and, if so, which convergents provide the 
solutions? Answers to both these questions are to be found in the next reading 
passage. 


(i) 


(i) 


(iii) 


READ B from the top of page 335 to page 337, line 3. 


Page 335 line 10 

The Lemma tells us which convergents can provide the solutions of Pell's 
equation. They are the convergents which immediately precede the 

partial denominator 2a, at the end of the period of the continued fraction 


for Ja: 


Jd = [403 41, 02,..., a2, d1, 2d | 


= [403 41, 45, . .., 42, 41, 2d, 41, A2, ..., A1, 20, 01, . .., 41, 2do,...]. 


The convergents corresponding to the arrowed partial denominators are 


ге, k = 1,2,3,..., where n is the length of the period. If n is even, the 
kn—1 


Lemma will tell us that pz, — d q2, ; = 1, giving a solution of Pell's 
equation for every positive integer k. On the other hand, if п is odd the 


Lemma gives pz, 1 — d qz, , = (—1), so that we have solutions for 
k = 2,4,6,... ; whilst for each odd value of k, the convergent gives a 
solution of x? — dy? = —1. 


Page 336, line —8 

Remember that the list above starts with Co = po/qo, and so p3/q; is the 
fourth in the given list. 

Page 337, line 3 

In Example 13-8, as the period of /13 has odd length, 5, it is only every 


tenth convergent, Pe. Pao Pao; ..., Which satisfies p? — 1342 = 1. However, 
do dio 429 
the Lemma on page 335 tells us that every tenth convergent, P» Bik 
44 dı4 
Еи satisfies p? — 1342 = —1. We can confirm that this is so for the 
24 


fourth convergent 18: 182 — 13-5? = —1, A 


B page 341, Problems 13.5, number 5. 
B page 341, Problems 13.5, number 4. 
B page 341, Problems 13.5, number 9. 


(a) 
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JB = [4;1,3,1, 8] 
As the length of the period is 4, which is even, solutions of Pell's equation 
Pak-1 


are given by convergents ‚ for all к > 1. The solutions we require will 
dak-1 


come from p3/q3 and p;/q;. The convergents are 


4 5 19 24 211 235 916 1151 
11 
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So solutions are (24,5) and (1151,240). 
(Check: 24? — 23-5? = 1 and 1151? — 23-240? = 1.) 


(b) /26 = [5;10] 


The period has odd length, 1, so the solutions to Pell's equation are given by the 


convergents Рок-1 for all k > 1. The first two solutions come from pi 3L 
42к-1 qı 10 
1 
(check: 51? — 26-10? = 1), and A = R (check: 5201? — 26-1020? = 1). 
3 


(c) 4/33 = [5;1,2,1, 10] 
The period length is 4 (even), so the first two solutions are рз/аз and p;/q;. The 
convergents of \/33 arej,$, 17, 22 211 270, 787 1057 . . So solutions are (23,4) 
and (1057,184). 
(Check: 23? — 33-4? = 1 and 10572 — 33-1842 = 1) в 

Suppose n + 1 = x? and + 1 = y*. Eliminating n from this pair of 


simultaneous equations gives x? — 2y? — — 1. Conversely, if x? — 2y? = — 1, 
then x? is odd and therefore can be written as x? =n + 1, where n is even. 
5 + 1. Since Jå = [1;2] has period of odd length, the 
Lemma on B page 335 guarantees infinitely many solutions (arising from the 
even convergents). The first solutions are (1,1) giving n = 0, (7,5) giving n = 48 
and (41,29) giving п = 1680. BN 


Solving for y gives y? = 


If ра then x? — dy” = —1 implies that x? = —1 (mod p). But —1 is not a 
quadratic residue of p = 3 (mod 4), and so no such x can exist. W 


On pages 337-9, B sets up a method by which all solutions of Pell’s equation 
can be generated from the least positive solution. We leave this as optional 
reading. However, to conclude the section, 


READ B page 340. 


The Pythagorean Equation 


You will be familiar with the fact that there is a right-angled triangle with sides 
of lengths 3, 4 and 5 (inches, say), which in terms of the celebrated Pythagoras’ 
Theorem says that 3? + 4? = 52. This provides us with a first solution of the 
Diophantine equation 


x-y-z (1) 
In this section we are going to find all solutions of this equation. We shall 
overlook the trivial solutions in which x or y is equal to zero, and, as x? — 


x? = (— х)? etc, we shall concentrate on solutions in which x, y and z are 
positive, since all other solutions are given by simple changes of sign. 


The forementioned solution leads immediately to infinitely many solutions of 
(1), since for any integer k > 1, 


(3k)? + (4k)? = (5K). 


Yet, in a sense, all these solutions are ‘essentially the same’—geometrically, the 
corresponding right-angled triangles are all similar. For this reason we attack 
the Diophantine equation (1) by first looking for positive solutions in which 
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ged (x,y,z) = 1. If we can find all such solutions, then all remaining positive 
solutions are found, as illustrated above, by multiplying each of the x, y and z 
values by any integer k > 1. 


READ B Section 11.1, page 242 as far as page 247, line —6. 


(i) Page 243, line 16 
As we can confine ourselves reasonably to the positive solutions of the 
Pythagorean equation, we shall also insist that the three integers in a 
primitive Pythagorean triple are positive. 


(ii) Page 244, line 13 
In any primitive triple, exactly one of x and y is even. Notice that if 
x? + y? = 22 then y? + x? = z? so that the values of x and y are 
interchangeable. To overcome this ambiguity, we adopt the convention that 
x takes the even value in any primitive triple. Note this convention carefully 
since it will be employed many times in what follows. 


(iii) Page 244, line — 13 
Having seen that in any primitive triple x is even, for it to be prime requires 
that x — 2. But we can easily show that there is no such triple, for 


22 фу? = 2 
gives 4 = (z — y) (z + y). 


As y and z are known to be odd, z — y and 2 + y are both even positive 
integers, giving z — y — 2 and z 4- y — 2. That is, z = 2, у= 0—a 
contradiction. 

On the next line, B is technically correct but has overlooked his convention 
that x is always even. He should say: ‘there are primitive triples in which y 
and z are prime’. 

There is overwhelming numerical evidence supporting the conjecture that 
there are infinitely many such triples, but, like so many similar situations 
that we have encountered in the course, this still awaits a general proof. 


(iv) Page 245, line 5 
The proof of Lemma 2 is completed more easily by using the result of 
SAQ 4 of Unit 2, which says that an integer is an nth power if and only if 
each exponent in its prime factorization is a multiple of n. Since ab is an nth 
power each of its exponents k,,..., K, Ji. j, is a multiple of n, and the 
expressions for a and b at the foot of page 244 confirm that they are indeed 
nth powers. 


(v) Page 246, line —5 
The result of Theorem 11-1 amounts to a complete classification of the 


primitive Pythagorean triples with x even. To get all positive solutions of 
the Pythagorean equation we drop the condition gcd (x,y,z) = 1, and have 


x = 2kst, y-—k($?—1?, z= k(s? + 22) 


for any integer К > 1, observing that x and у are interchangeable. A 


B page 249, Problems 11.1, number 1. 


Show that in any solution of x? + y? = z? at least one of x, y or zis a multiple of 5. 
Hence solve: 


B page 249, Problems 11.1, number 4. 
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SAQ 7 Show that if p is an odd prime, then x? + y? = 22 has exactly one positive solution 
with y = p. 


SAQ 8 (i) Show that if N, a and b are positive integers such thàt N — ab, then, for 
any integers s and t, 


х= as” — bt’, y 22st, z= as? + bt? 
is a solution of x? + Ny? = z2, 
(ii) There are six positive solutions of x? + бу? = z? in which z < 20: find them. 


SOLUTION (а) Theorem 11-1 tells us that primitive solutions are given by 


TO BARES х= 25, у= 52-1, z=% 42? 


where gcd (5,1) = 1, s > t > 0 and sz t (mod 2). Given that x = 16, the only 
primitive solution arises from s = 8, t = 1, whence x = 16, y = 63 and 
2 = 65. 


We can find solutions which аге not primitive by the same method but 
solving for x being a divisor of 16. 


s=4,t=1 gives х= 8, у = 15, z = 17. 

Multiplying by 2 gives the solution x = 16, у = 30, 2 = 34. 
s=2,t=1 gives х= 4, у= 3,2= 5. 

Multiplying by 4 gives the solution x = 16, у = 12, 2 = 20. 


(b) For x = 40 we require all pairs of positive values of s and t (s > t) such that 
2st = 40, gcd (s,t) = 1, s # t (mod 2). They are 


s=20,t=1 giving x = 40, y = 399, z = 401, 
and 
s=5,t=4 giving х= 40, у= 9,2 = 41. 
For x = 60, proceding similarly, 
s=30,t=1 giving x 60, y = 899, z = 901, 
s=15,t=2 giving x = 60, у = 221, z = 229, 
s=10,t=3 giving x=60,y=91, z= 109, 
$—6, t=5 giving x=60,y=11, z=61. B 


SOLUTION If x? + y? = 22, suppose that 5x and 5y. We must show that 5lz. 
Eo SAQ 6 Since 5x, x? = 1 or 4 (mod 5) and similarly y? = 1 or 4 (mod 5). Therefore 
x? + y? = (1 or 4} + (1 or 4} (mod 5) 
= 0, 2 or 3 (mod 5). 


But 2 and 3 are not quadratic residues of 5 and so we must have z? = 0 (mod 5), 
which amounts to 5|z. 


To prove that 60|xyz, we show that 3|xyz, 4|xyz and 5|xyz, and appeal to 
Corollary 2 on B page 28, observing that 3, 4 and 5 are pairwise relatively 
prime with product equal to 60. In fact, we have just shown that 5|xyz and on 
page 247 of B we saw that 3|x or 3]у, yielding 3|xyz. So it remains to show only 
that 4|xyz. For this we need only look at the form of x in any primitive 
solution. By Theorem 11-1, x = 2st where s Æ t (mod 2). The latter fact forces 
one of s or t to be even, and in turn 4|x. Therefore 4|xyz, as required. Bl 
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Suppose there is a solution with y = p. Then from x? + p? = z? we have 
p = 2 - х? = (z х)(2 + х). 


Now x and z are both positive and z > x, so z + x > z — x > 0. As p? can be 
expressed as a product of two positive integers in only two ways, namely 1-р? 
апа p-p we conclude that 

z—-x=landz+x=p’. 
(That both should be equal to p is excluded since z + x > z — x.) From these 
two equations we conclude that 


lord 1.2 
Ы = + 1). 
x zP 1) and z zP ) 


So any solution with у = р is uniquely determined. On the other hand, it is a 
simple matter to verify that 


1 
1) y= ,م‎ 2 = 5(р? +1) 
is a primitive Pythagorean triple, for any odd prime р. 


COMMENTS 


1. In fact, it is not difficult to show (from Theorem 11-1) that not only odd 
primes, but every integer п > 1 (provided n # 2 (mod 4)) occurs in some 
primitive Pythagorean triple. The reason why п = 2 (mod 4) is excluded was 
seen in SAQ 6, where we found that the even member in any primitive 
Pythagorean triple is always divisible by 4. 


2. If z= $(p? + 1) is prime then we have found a primitive Pythagorean triple 
with y and 2 both prime. We commented in Note (iii) that it is not known 
whether there are infinitely many such triples. Our solution therefore brings 
to the fore yet another unsolved problem: are there infinitely many pairs of 
primes p and q such that 2q = p? +17 Ш. 


(i) This part is just algebraic manipulation: 
x? + Ny? = (as? — bt?)? + ab(2st)? 
= a’s* + 2abs?t? + b?t4 
= (as? + bt?)? 


= z?, 


(ii) For N = 6 we have the two cases a = 6, b = 1 and a 3,b=2 to 
consider. (The cases a = 1, b = 6 and a = 2, b = 3 will lead to the same 
solutions but with a change of sign in the x-value.) 


With a = 6 and b = 1, z = 6s? + t? < 20 requires s = 0 or 1. However, 
s = 0 will give y = 2st = 0, so we are left with s = 1, in which case t = 1,2 
or 3: 


$=] tsi KS ya 2. жш, 
s=Lt=2 => х=2,у= 4,2 = 10, 
$—1t1-23 = х= –-3,у= 6,2 = 15. 


As the sign of x is immaterial, since it is to be squared, the last gives us the 
positive solution x = 3, у = 6, = 15. 


With a = 3 and b = 2, z = 3s? + 212 < 20 gives us three cases: 


76 


8.3 


NOTES 


SAQ 9 
SAQ 10 


M335 8.2/8.3 


s=l,t=1 => x=ly=2,7=5, 


s=1t=2 => х=—5,у=4,л= 11 


(this gives the positive solution x = 5, у = 4, z = 11), 


Fermat’s Last Theorem 


In this section we shall read more about the history and the present state of 
play of attempts to prove Fermat’s Last Theorem. We shall also see a proof for 
the particular case n = 4. In that proof, we shall encounter an important and 
powerful technique known as ‘Fermat’s method of infinite descent’. It is used 
here to show the non-existence of solutions of particular Diophantine equations. 
In Section 8.4, however, we shall use the descent method in a slightly different 
way to establish the existence of solutions of certain other Diophantine 
equations. 


As we shall use it in this section, the idea behind the descent method is as 
follows. We suppose that the Diophantine equation in question has a positive 
solution. Then, focussing on one of the variables involved, say z, the Well 
Ordering Principle tells us that amongst all solutions there is one in which z 
takes its least positive value. We contradict this by showing, by some algebraic 
argument, that any solution of the Diophantine equation gives rise to a second 
solution in which the positive value of z has decreased. (This key stage of the 
method is called the descent step.) The only conclusion to be drawn is that the 
set to which we have applied the Well Ordering Principle must be empty. That 
is, the Diophantine equation has no positive solutions. 


The application of the descent method in the reading passage is quite involved, 
but in the SAQs at the end of the section we have given you two simpler 
Diophantine problems which can be solved by this method. 


READ B Section 11.2, page 250 as far as page 254, line 16. 


(i) Page 251, line 4 
Since xo = d x, and yo = dy, we have xé + y$ = d* (xt + yt) = 22. From 
this it follows that d?|zo, which justifies putting 20 = d?z,. 

(ii) Page 251, line 19 
1 = y$ (mod 4) follows since yo is presumed to be odd and the square of 
any odd integer is congruent to 1 modulo 4. 

(ili) Page 254, line 16 
It is said that the great number theorist Landau had postcards printed 
which read: ‘Dear Sir or Madam: Your attempted proof of Fermat’s 
Theorem has been received and is returned herewith. The first mistake is 
оп page line . Landau left to his students the task of filling in the 
missing numbers. A 


B page 258, Problems 11.2, number 12. 
Use the method of infinite descent to show that 
х? + 3ر3‎ = 923 


has no positive solutions. (Hint: In any solution, 3|x.) 
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Show that if x^ + y? + z = 0 (mod 4), then each of x, y and z is even. Hence 
use the method of infinite descent to show that 


x? +y? +22 = 2xyz 
has no solutions other than x = y = z = 0. 


[Hard] B page 258, Problems 11.2, number 7. (Hint: The hint in B provides 
the final step of the argument but it is far from clear how one arrives at his 
expressions for х? + y?, x + y and x — y. Observing that the product of these 
three terms is x^ — у“, Lemma 2 (B page 244) will give us the expressions, 
provided the three factors are pairwise relatively prime. For that we first need 
gcd (x, y) = 1.) 

NT 
x + ж = 5 rearranges to x* + y* = Е 2 
this equation has по solutions in positive integers. IN 


2 
| ‚ and Theorem 11-3 tells us that 


Suppose we have a positive solution Xo» Jo, Zo SO that 
XQ + Зуў = 922. 

Now хў = 3(328 — yê), showing that 3|xp. Putting хо = 3x,, 

27x1 + 3yê = 92%, or 9x? +y? = 322. 
But now 3|y, and so, putting yo = 3y,, 

9х? + 27у? = 325, or 3x} + 9y3 = 25, 
This time 3[zo, and putting zo = 3z, gives 

3x1 + 9y? = 2723, or x} + 3у = 923. 


This completes the descent step. Any positive solution gives rise to a another 
solution in which the value of each of the variables has decreased, for example 

x "m А : 
О& = F < xo. Hence, by the method of infinite descent, the Diophantine 
equation has no positive solutions. 


COMMENT In this example we could avoid appealing to the method of 
infinite descent. Assuming the existence of a solution, cancellation of all 
common factors would give a solution with ged (x,y,z) = 1. The above 
argument produces a contradiction by showing that gcd (x,y,z) > 3. Ml 


а? = 0 or 1 (mod 4) according as a is even or odd respectively. So 
x? + y? + 2? = (0 or 1) + (0 or 1) + (0 or 1) (mod 4) 
= 0, 1, 2 ог 3 (mod 4) 


according as none, one, two, or three, respectively, of x, y and z are odd. So 
2 2 v — я, 
X^ + y^ + 2° = 0 (mod 4) precisely when x, y and z are all even. 


In looking for solutions of x? + y? + z? = 2xyz, we assume that x, y and z are 
positive, since all other non-trivial solutions are given by simple changes of sign. 


If x, y and z are all odd then x? + y? + z? = 3 (mod 4) whilst 2xyz = 2 (mod 4). 
So at least one of x, y and z is even, and now 2xyz = 0 (mod 4). By the first part 
of the question, it follows that x, y and z must all be even. Therefore, writing 

x = 2х', y = 2y' and z = 2z', we have 


Qx'Y + Qy'Y + Qz'yY = 20x?) (2y’) (22’) 


x2 + у + 2ے‎ = 4Ax'y'z. 


As the right-hand side is still congruent to 0 modulo 4, x', y' and z' are all even 
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(this by the first part again). Writing x’ = 2х”, etc., we now obtain 


xu y" + = 8х”у”т” 


and yet again х”, y" and z” are all even. We could continue forever in this way, 
halving the x, y and z values and yet retaining even, positive integers; for from 


xê + уф + 28 = 2'xoyozo (r z 1) 
we have 
2 2 2 
Xo Уо 20 r*1[Xo| [Уо | [Zo 
+ + = 
aq. Xo" Y6 Zo PN 
with rues and > even. Hence, by the method of infinite descent, the 


Diophantine equation has no solutions. Hi 


In looking for solutions of x* — y* = 222, we may assume, without loss of 
generality, that gcd (x, y) = 1; for, if gcd (X, y) = d > 1, then putting x = d x,, 
y — dy, and z — d?z, (since the equation shows that d?|z) gives xf — y = 222, 
where рса (х,у) = 1. 


As x* — y* is even and ged (x, y) = 1 it follows that x and y are both odd, and 
in consequence that x? + y?, x + y and x — y are all even. So x? + y? = 2r, 
x + y = 2s and x — у = 21, for some integers r, s and t. Now 


xt — yt = (x? + y?) (x + у) (х y) = 8rst = 2z?, 


Z\2 
d t= {-]. 
апа SO rs: B 


Suppose for the moment that r, s and t are pairwise relatively prime. Then 
Lemma 2 (B page 244) tells us that r, s and t are each squares, so that 

x? + y? = 2a?, x + y = 2b? and x — y = 2c? for some positive integers a, b and 
c. From this we have 


4b* + 4c* = (x + у)? + (x — у)? = 2x? + 2y? = 442, 


But the equation bt + c^ = a? is known to have no positive solutions (Theorem 
11-3), contradicting the existence of such integers x and y. 


To tidy up the loose ends we must finally confirm that 
gcd (s, t) = ged (r,s) = ged (r,t) = 1. 


If dls and dlr then d|s + t. That is, d|x and ау. But ged (x,y) = 1, and so 41. 
Therefore gcd (s, t) = 1. 


If pir and p|t (where p is prime), then p|r — 212. That is, р|ху and so by Euclid's 
Lemma ріх or ply. But p|x — y (=2t), and so р|х and ply, which contradicts 
ged (x,y) = 1. 


Similarly, if pr and p|s then p|2s? — r. That is plxy, which as above leads to a 
contradiction of рса (x, y) = 1. IN 


Sums of Squares 


Chapter 12 of B is devoted to the problem of representing integers as sums of 
squares. The culmination of the chapter is a theorem, due to Lagrange, which 
proves that every positive integer can be expressed as a sum of squares of four 
integers; in other words, the Diophantine equation 


N=w? +x? +y +22 


has solutions for every positive integer N. 
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Since we allow the variable to take value zero, we have, as an immediate 
Corollary of Lagrange’s Theorem, that every positive integer can be expressed 
as a sum of n squares, for each n = 4. So, in addition to the above result, there 
are just two interesting cases: which integers can be expressed as a sum of two 
squares, and which as a sum of three squares? In this section we shall give most 
of our attention to the first of these problems, which amounts to finding those 
integers N for which the Diophantine equation 


N=x? +y 


is solvable. Thereafter we shall sketch a proof of Lagrange’s Theorem, and take 
a brief look at the three squares problem. 


READ B Chapter 12, page 260 as far as page 263, line 8. 


Use the Lemma to express 4420 as a sum of two squares. 


4420 = 22:5 -13 -17 = (2? + 02) (2 + 12) (3? + 22) (4? + 12). 
Using the formula on the final line of page 262 of B three times: 
4420 = (4? + 27) (3? + 27) (4? + 12) 
= (16? + 2?) (4? + 12) 
= 667 + 82. 
COMMENT The answer is not unique. Indeed, in addition to the above, 
4420 = 64? + 18? = 62? + 24? = 48? + 462. 
We obtain the other solutions by changing the order of the bracketed terms, or 


by changing the order of the two squares within the brackets. For example, 
from the final line of the above solution, 


4420 = (16? + 27) (12 + 42) = 24? + 622. Ш 


To investigate whether or not the integer N is a sum of two squares, the Lemma 
suggests that we should look at its prime factorization. If each prime divisor of 
N is a sum of two squares then so too is N itself. Now primes come in 3 types: 


(i) the even prime, 2 = 1? + 17, is a sum of two squares; 

(ii) odd primes of the form 4k + 3, which, as we have seen in Theorem 12-1, 
cannot be expressed as a sum of two squares; 

(ii) odd primes of the form 4k + 1. 

We shall now prove that all primes of type (iii) are sums of two squares, and 

having done this it will be a simple matter to complete the two squares 

problem. We are going to give an alternative proof to the one in B, our proof 

being essentially B's proof of Lagrange's Theorem—the details are different, of 

course, but the methods are identical. Once you have sorted out the proof given 


below you should have no difficulty with B's proof of Lagrange's Theorem 
(which we have left as optional reading). 


The proof uses Fermat’s descent method, but in a different way from our 
previous applications. 


(Fermat) 


If p = 1 (mod 4) then p can be represented as a sum of two squares. 
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Preamble to the proof 


We start by assuming that some multiple of p can be represented as a sum of 
two squares. More precisely, we assume that 


mp =x? + y? (1) 
has a solution for some 1 « m < p. From this we deduce that a smaller multiple 
of p can be so represented: 

m'p = u? +1? 
where 1 x m' « m. This is our customary descent step. 


The difference is that this time we shall show that (1) does have a solution. The 
descent step tells us that there is then a smaller solution, and from this a smaller one, 
and so on. The process cannot go on for ever, and the only way it can terminate is by 
descending to a solution with m' — 1 (so that the descent step cannot be applied 
again). The inescapable conclusion is that p = x? + y? has a solution. 


And now for the details. 


Part (i) (The descent step) 
Suppose that the equation mp = x? + y? has a solution in which 1 < т< p. Let 
u and v be the least absolute residues modulo m of x and y, respectively. That is, 


u=x(modm), v=y(modm); -—im<uv< im. 
Then 
u? v? = x? + y? = 0 (mod m), 
and so 
u? + vu? = mr (for some r > 0). 
Now if r = 0, then u = v = 0 implying that m|x and mly. But then, from 
mp = x? + y?, we conclude that m|p—which is a contradiction. So, 


1 lím m? 
O<r=—{u? +97} < —{— 4 —_ ; 
r mit est urbi 
Multiplying together the equations mp = x? + y? and mr = и? + v? gives 
mrp = (x? + y?) (u? + v?) = (хи + yo)? + (xv — yu), 
Now 
хи + yo = х? + y? =O(modm), implying that т|хи + yo, 
and 
Xv — yu = xy —xy=O0(modm), implying that m|xo — yu. 


Putting xu + yv = mX and xv — yu = mY gives 
mrp = т?Х? + mY? 
or rp = X? + Y?, where 1 < r < m. 


Part (ii) (To show mp = x? + y? has a solution with 1 « m <p) 
Since —1 is a quadratic residue of p (Corollary, B page 193), the equation 
x? + 1 =0 (mod р) has a solution г. Selecting the least positive residue, we may 
assume that 0 < z < p — 1. So there is an integer m > 1 such that 
mp = 22 + 1? 
and 


1 1 
m= {2 +17} 5210 - 1y + 1) Pis = 2(р – 1) <p. 


So there is a solution (1 < т < р) to mp = х? + y? and hence, by part (i), а 
solution to х? + y? = р for p = 1 (mod4). W 
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Given that 60? + 1? = 13-277, retrace the proof of the descent step to express 
the prime 277 as a sum of two squares. 


From 

60? + 1? = 13-277, (1) 
we have т = 13, x = 60 and у = 1. So и = 60 (mod 13) and v = 1 (mod 13) with 
—7 <u, v <7, giving u = —5 and v = 1. Therefore 

(—5 +1? = 2-13 (2) 


so that r = 2. Multiplying equations (1) and (2) together, 
2-13?-277 = (60? + 12) ((—5)? + 12) 
= (—299)? + 652, 
and on dividing both sides by 13? (and discarding the minus sign inside the 


square), 
sin 2-271 223? + 52, 


completing the first descent step. 
Repeating with m — 2, x — 23 and y — 5 we have u — v — 1 and from 
и? +? = 1? 1? — 1:2, 
r = 1. Multiplying the two equations, 
27۰277 = (23? + 5?) (1? + 17) = 28? + 182, 
and on dividing both sides by 22, 
277 = 14° +92. WI 


After working through SAQ 14 you may feel that there ought to be an easier 
way of expressing an integer (like 277) as a sum of two squares. Indeed, you 
would be justified in thinking that a ‘trial and error’ method would get you to 
the solution more easily. Knowing that 277 can be expressed as a sum of two 
squares, we need only check 277 — 12, 277 — 22,..., etc. until we find a square, 
in the knowledge that one must turn up by the time we reach 277 — 122 (since 
12? > 234), 


Knowing that a prime of the form р = 4k + 1 can be represented as х2 + у?, 
mathematicians understandably turned to the task of finding ways of 
constructing the integers x and y in terms of p. Several such constructions 
emerged, but none of them particularly easy to employ. The first, due to 
Legendre, showed how to obtain x and y from the continued fraction of Jp. 
The disadvantage of this method is that we first have to find the continued 
fraction of JP, which is no mean task when p is large. The next, due to Gauss, 
is easy to state: 

2k) 


! 
If p = 4k + 1, then x = Жүр mod р) and у = (2k)! x (mod p), 


where x and y are least absolute residues modulo p, satisfy x? 4- y? — p. 


This is very neat, but try using it to express 277 (k — 69) as a sum of two 
squares! 


For practical purposes, undoubtedly the best method for representing an integer 
as a sum of two squares is to resort to the suggested trial and error method, 
coupling it, of course, with the construction given in the Lemma at the foot of B 
page 262. 
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Express 5321 = 17-313 as a sum of two squares, 
Testing 313 — 17, 313 — Q^ ees (OF being a square, we eventually find that 
313 — 12? = 132. Therefore 313 = 122 + 13? and 
5321 = (4? + 12) (122 + 132) = 61? + 402. 


Had we written 5321 = (4? + 12) (13? + 12°), the formula for this product at the 
foot of B page 262 gives the alternative answer 


5321 = 64? + 352, 
These two are the only solutions in non-negative integers. Ш 
We have seen that some integers can be represented in more than опе way as a 


sum of two squares. This is not true of primes. A prime of the form 4k + 1 is 
representable in a unique way as a sum of two squares of positive integers. 


READ B page 265, from line 7 to the foot of the page. 


Page 265, line 12 
The conclusion of this line is justified by 
ad? — b?c? = (a? + b?)d? — (e? + d^)? = pd? - pb. A 


Having solved the problem of which primes are sums of two squares we are 
now ready to determine which integers are sums of two squares. 


READ B from page 267, line 7 to page 269, line 15. 


Which years in this decade (1980-1989) can be represented as a sum of two 
squares? For each affirmative answer find one such representation. 


B page 271, Problems 12.2, number 3 parts (a), (b), (c) and (d). 


Show that if p = 4k + 1 is prime, then 2p = x? + y? has exactly one solution 
with 0 < y < x. 


For which integers n(0 < n < 35) is the congruence 
x? + y? = n (mod 36) 


not solvable? Hence show that at least 40 per cent of the positive integers < 10% 
are not expressible as a sum of two squares. 


1980 = 2?-3?-5-11 1985 = 5-397 
1981 = 7-283 1986 = 1 
1982 = 2-991 1987 = prime 
1983 = 3-661 1988 = 2?-7-71 
1984 = 2°-31 1989 = 32-13-17 


Excluding 1985 and 1989, each of the remaining years is divisible by a prime 
р = 4k + 3 but not divisible by p?. By Theorem 12-3, these years are not 
representable as sums of two squares. 


1985 = 32? + 31? = 442 +72, 
1989 = 332 + 302 = 422 4152, gm 


SOLUTION 
TO SAQ 17 


SOLUTION 
TO SAQ 18 


SOLUTION 
TO SAQ 19 


(a) Since 2" is not divisible by any odd prime, the resuit follows from Theorem 
12-3. More specifically: 
n even => 2" = Q"?y + 02, 
n odd => 2" = 27272 72, 
(b) Modulo 9, the squares are 0, 1, 4 or 7. So for any two integers x and y, 
x? + y? = (0, 1, 4, 7} + {0, 1, 4, 7) 
= 0, 1, 2, 4, 5, 7 or 8(mod9). 
Therefore n = 3,6 (mod 9) cannot be represented as the sum of two squares. 


_ r(r-- 1) , s(s +1) 
UM. v р 


4n +1 —2r(r +1) + 2s(s +1) +1 = (r +s +1) + (r— sy. 
(d) F, = 22") + 12, for n >1. Ш 


then 


(c) Ifn 


If 2p = x? + y?, then x? + y? = 2 (mod 4), which means that x and y are both 
odd. In that case 


+y\? [x —y)\? 
5 02901 
х +у 2 + 2 r 
х+у\2 [x —y\? 
so that p = 2 T xj But we know that p can be represented as a 
sum of two (positive) squares in a unique way, and so i = and * 7 у (апа їп 


consequence x and y) are uniquely determined. Ё 


x? + y? = п (mod 36) has solutions if and only if each of x? + y? = n (mod 4) 
and x? + y? = n (mod 9) have solutions. We saw in SAQ 17(b) that 

x? + y? =n (mod 9) has solutions except for n = 3 and 6. Similarly, x? + y? =n 
(mod 4) has solutions for n — 0, 1 and 2 but not for n — 3. It follows that 

x? + y? =n (mod 36) has solutions for all but the following fifteen values of n: 


п = 3, 6, 7, 11, 12, 15, 19, 21, 23, 24, 27, 30, 31, 33, and 35. 


Consequently, in each set of 36 consecutive integers there are at least 15 which 
are not sums of two squares. (Note that if n is congruent modulo 36 to a sum of 
two squares, it does not follow that n itself is a sum of two squares. For 
example, 38 is not a sum of two squares, but 38 = 1? + 1? (mod 36).) Therefore 
the percentage of positive integers not exceeding 10° which are not sums of two 
squares is at least 


4-106 — 36 

36 Е m" 

TIU © 100 = 41.7 > 40. 
(Note that since 36 does not divide 10° exactly, we subtract 36 in the numerator 
to play safe.) Bi 


To round off the section we shall take a brief look at the three squares and four 
squares problems. 


As commented earlier, the proof in B of Lagrange’s Theorem, that every positive 
integer is a sum of four squares, is remarkably similar to the proof we have 
given for the two squares problem. It hinges around an identity (Lemma 1, B 
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SOLUTION 
TO SAQ 20 


M335 8.4 


page 275) showing that if two integers are each representable as a sum of four 
squares then so too is their product. This suggests that we attack the problem 
by establishing that every prime is a sum of four squares. The even prime 

2 = 17 + 1? + 0? + 0? creates no problem. The proof for the odd primes uses 
the descent method in just the same way as in the case of two squares. 


For any odd prime p, we assume that 
mp = м? +x? + y? + 2? 


has a solution for some m such that 1 « m « p. We achieve the descent step by 
deducing from this that 


m'p=s? +t? + и? +1? 


for some smaller multiple m' such that 1 < m' < m. The second part of the proof 
shows that the initial assumption, that an appropriate multiple of p can be 
expressed as a sum of four squares, is valid. (This is Lemma 2, B page 275.) 


The three squares problem is much more difficult, simply because it is not true 
that if two integers are each the sum of three squares then so too is their 
product. For instance, 3 = 1? + 1? + 1? and 5 = 2? + 1? 4 02, but their 
product, 15, is not a sum of three squares. This forces a completely different 
approach to the problem. 


Discovering which integers cannot be represented as a sum of three Squares is 
an easy matter; the difficulty comes when one tries to prove that these are the 
only such integers. 


READ B Section 12.3, page 273 as far as page 274, line —3. 


Determine whether or not each of the following integers can be represented as a 
sum of two squares and/or three squares. 


(a) 1,000,003 (prime), (b) 10!, (c) 100! (Hint: Use Theorem 6-9, B page 126). 


(a) 1,000,003 = 3 (mod 4), so it is not a sum of two squares. 
1,000,003 = 3 (mod 8), so it is a sum of three squares. 


(b) 10! = 1:2:3-...- 10 = 25 3*- 5? - 7. Since 7]10! but 72/101, 10! is not a sum 
of two squares. 
Writing 10! = 4* - (3* -5?- 7), and discovering that 3*-5?-7 = 7 (mod 8) 
(which is easily seen since 3? = 5? = 1 (mod 8)), we have that 10! is of the 
form 4* (8m + 7) and so, by Theorem 12-5, is not a sum of three squares. 


c) Since 100! is divisible by several primes p = 4k + 3 but not by their squares 
y 
(e.g. p — 59), it follows that 100! is not a sum of two squares. 


Using Theorem 6-9 (B page 126), the exponent of the highest power of 2 
dividing 100! is 


100 á 100 ái 100 " 100] [100 " 100 
7 4 8 16| |32 64 
= 50 +25 +12 +6 +3 +1 = 97. 


Therefore 100! = 2?" т = 4*8 (2m), where m is odd. As this is not of the 
form 4" (8m + 7), 100! is a sum of three squares. В 
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Ѕиттагу 


The study of Diophantine problems is a vast subject lying at the heart of 
number theory. In this unit we have looked at a selection of historically 
significant Diophantine equations. 


Section 8.1 investigated Pell’s Equation, x? — dy? = 1, and the closely related 
equation x? — d y? = —1, where d > 2 is a non-square integer. We discovered 
that all solutions of these equations are readily obtained from the convergents of 
the continued fraction of Ja. 


In Section 8.2, we focussed on the equation associated with Pythagoras’ 
Theorem, x? + y? = z?. The geometric interpretation of the corresponding 
Diophantine problem is to find all right-angled triangles with sides of integer 
length. Appreciating that all positive solutions arise as scalar multiples of 
primitive ones, i.e. ones in which gcd (x, y,z) = 1, we successfully classified all 
primitive Pythagorean triples. 


Section 8.3 centred around Fermat’s Last Theorem, the unresolved conjecture 
that x" + у" = 2" has no non-trivial solutions for п 2 3. Using our result on the 
Pythagorean equation, we proved Fermat’s Last Theorem for the particular case 
п = 4. The strategy of the proof incorporated Fermat's method of infinite 
descent, for which we subsequently found several other applications. 


Section 8.4 turned to the problem of representing integers as sums of squares. 
We concentrated on the two squares problem, classifying those integers N for 
which the Diophantine equation N — x? 4- y? has solutions. We observed that 
Lagrange's Theorem, which states that every positive integer is a sum of four 
squares, can be proved in a very similar manner to the two squares problem, 
and we outlined that proof. Finally, we discovered that no integer of the form 
4* (8m + 7) can be expressed as a sum of three squares. We omitted the 
(difficult) proof that all other positive integers are sums of three squares. 


The main definitions and results encountered in this unit are listed in the 
corresponding entry in the Handbook. 
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